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I. INTRODUCTION 



The Standard Model of particle physics (SM) contains six quarks that can be arranged into three generations of 
two quarks each, and six leptons, similarly arranged. The top quark is the wcak-isospin partner to the bottom quark, 
which together constitute the third generation of quarks. The existence of three quark generations was postulated 
as early as 1973 by Kobayashi and Maskawa since mixing am ong three generations (and no fewer) could provide 
a mechanism for CP violation (IKobavashi and Maskawal 119731 ). E xperimentally t he multi-generation quark m odel 
found firm footing in 1974 with the discovery of the charm quark ( Aubert et al. , 1974 lAugustin et all 1974 ). the 
weak-isospin partne r to the strange q uark in the second generation of quarks. The discovery of the third generation 



tau lepton in 1975 (IPerl et 



discovered in 1977 (Herb et 



1975) implied the existence of a third generation of quarks. The bottom quark was 



19771 ) and the search for its weak-i sospin partner began in earnest . It was nea rly two 



decades before the top quark was directly observed by experiment ( Abachi et al. , 1995t Abe et al. , Il994 1"l995l ) owing 
to its surprisingly large mass. Of course the immediate consequence of this discovery was to further solidify the SM. 
Indeed, by the time the dire ct observation was made th ere was mounting indirect evidence that the top quark existed 
and that it had to be heavy ( Schaile and Zerwasl 1992 ). But it could be argued that this indirect evidence was model 
dependent and that an important consequence of the top-quark discovery was the definitive elimination of "top-less" 
theories that might otherwise still be viable. In fact, the principal interest in the top quark arises because it offers a 
window into new physics, into physics Beyond the Standard Model (BSM). Thanks to its very large mass (it weighs 
about as much as a gold atom) BSM contributions to physics involving top quarks can occur in a wide variety of ways 
- the production mechanisms can be affected, the decay widths can be altered, its intrinsic properties changed, or the 
experimental signature mimicked by a new particle of similar mass. A thorough exploration of all these possibilities is 
what constitutes the field of "Top Quark Physics" , pioneered over the last 15 years by the experiments at the Fermilab 

Tevatron. _ ^ ^ 

I ' recent years there have been several reviews of top-quark phys ics ( Demina and ThomsonL 2008 ; Incandela et al. 
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20091 : iKehoe et al\ . l2008t IPleierl . 120091 : iQuadtl . 120061: IWagnerl . l2010l : |WickeJ, |2010|) . In this article we will concentrate 



on what has been experimentally established by the Tevatron experiments and the analysis techniques required to do 
so. The emphasis of this article will be on t he pp — > ti process responsible for most of the progr ess to date. While the 
pp —> tq process has been recently observed ( Aaltonen et al. , 2009d . 2010J; Abazov et aZ.L[2009b ) the available samples 
are still rather small, which limits their sensitivity relative to what can be accom plished with the ti sample. Recent 
revie ws of the emerging field of top-quark studies using the tq events are available ( Heinsonl . 2010tlHeinson and Junkl . 
20111 ). The study of top quarks is expected to evolve over the next years using data from CERN's Large Hadron 
Collider (LHC). Where appropriate we comment on differences between the Tevatron studies and the LHC possibilities. 



II. STANDARD MODEL PRODUCTION AND DECAY 



We begin by introducing some concepts and phenomenology to which we will refer throughout this article. We 
will not summarize the Standard Mo del theory in any detail. This has already been done many times and interested 
readers should consult the references ( Kronfeld and Quiggl 20101 : Quiggl 2009 ) for in-depth discussions. 

Simply put, high energy physics endeavors to enumerate the elementary particles which constitute all known matter 
and to characterize the interactions which govern their behavior. The best guess to date is called the Standard Model of 
particle physics. The SM is a quantum gauge theory based on the SU(3) gauge symmetry of quantum chromodynamics 
(QCD) and the SU(2) x £7(1) gauge symmetry of quantum electrodynamics (QED). 

In the SM there are twelve elementary particles - six quarks and six leptons, each arranged into three generations 
as shown here 



t 
b 

v r 

T 



(1) 



The up- type quarks, up (u), charm (c), and top (t), have electric charge + § while the down-type quarks, down (d), 
strange (s), and bottom (6), have electric charge —3. The neutrinos, v ei v^, and v T are the electric-charge-neutral 
lepton partners to the electron (e), muon (/i) and tau (r) leptons, respectively. The e, /1, and r all carry electric 
charge —1. The measured masses of these particles are given in Tab. HI All these particles have anti-particles with 
the same mass and quantum numbers except for the opposite electric charge, the opposite parity, and being singlets 
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rather than doublets of the ST/ (2) group. Often anti-particles are represented with the same symbol as the particle 
but with a bar; for example a quark is represented as q and an anti-quark as q. An important subtlety worth noting 
is that the weak eigenstates of the down-type quarks are mixtures of the mass eigenstates. This mixing is represented 
by the CKM matrix, 

V u d Vus V u b 

V cd V cs Vcb I I 8 I . (2) 





V td V ts Vtb 

The interactions of the elementary particles are governed by three forces, the strong force of QCD, and the electro- 
magnetic and weak forces of QED. These forces are mediated by the exchange of bosons. The strong force is mediated 
by the exchange of gluons (g), the electromagnetic force by the exchange of photons (7), and the weak force by the 
exchange of the charged W bosons (W ) or the neutral Z bosons (Z°). The strengths of these interactions are charac- 
terized by coupling constants with these approximate relative strengths a s trong : ot em : a W cak ~ jq '■ j^q '■ 10 ooo ■ While 
the gluon and photon are massless, the W and Z b osons are massive and have been measured to be mw = 80.40GeV/c 2 



and mz = 91.188 GeV/c 2 (jNakamura et al . 2010). Figure [T] shows Feynman diagrams for an example of strong, elec- 



tromagnetic, and weak interactions. 

Not all of the elementary particles participate in all three interaction forces. The photons couple to the electric 
charge, so neutrinos, for example, do not participate in the electromagnetic interaction. The gluons couple to the 
color charge. Only quarks carry color charge and participate in the strong interaction. All the quarks and leptons 
can participate in the weak interaction. The gauge bosons can interact subject to these same constraints. Gluons 
carry color and can self-interact. Similarly W and Z bosons can interact with each other. The electrically charged 
W boson can couple to the photon. 

In the simplest form of the SM, the SU(3) x SU(2) x U(l) symmetries predict that all the elementary particles 
and the gauge bosons are massless, which is clearly discrepant with experimental observations. The easiest way to 
remedy this is to introduce a scalar field whose vacuum state breaks the SU{2) x U(l) electroweak symmetry. With 
a non-zero vacuum expectation value, v = 246 GeV, this field, called the Higgs field, imparts mass to the W and 
Z bosons via gauge interactions, and to the quarks and leptons via Yukawa interactions. The electroweak gauge 
interactions establish a well defined relationship between the Higgs vacuum expectation value and the masses of the 
W and Z bosons. The quark and lepton masses are linearly proportional to a Yukawa coupling, Y q> £, and the Higgs 
vacuum expectation value. From their measured masses we can infer that the top-quark Yukawa coupling is of order 
Y t ~ 1, while the others are significantly smaller ranging from Y e ~ 10 -6 for the electron to Yf, ~ 10~ 2 for the 6-quark. 

The physical manifestation of the Higgs field is a scalar boson, H°. The Higgs boson is the only SM particle that 
has not been experimentally observed. While the mass of the Higgs boson is not predicted by the SM, it can be 
inferred since quantum loops including the Higgs boson induce corrections to several experimental observables at a 
level comparable to the precision with which these observables are measured. A global fit to these so called "precision 
electroweak observables" pr edicts that the mass of the Higgs boson lies in the range of about 40 < run < 160GeV/c 2 
at the 95% confidence level ( Flacher et al. . 20091 : Nakamura et al . 2010l ). Direct searches for the Higgs boson at LEP 



and the Tevatron have ruled-out the masses tuh < 115 GeV/c 2 and 158 < m# < 175 GeV/c 2 at the 95% confidence 
level. It should be noted that all these constraints are valid only in the context of the SM. The electroweak symmetry 
may well be broken in some other way. 

A peculiarity of QCD worth noting is that quarks can only be observed as bound state hadrons. There are two types 
of hadrons: mesons are quark-anti-quark bound states, and baryons are three-quark bound states. An important 
experimental consequence of this is that quarks produced in a high energy pp interaction manifest themselves as 
collimated streams of hadrons called "jets". The energy and direction of a jet are correlated to the energy and 
direction of its parent quark. The process by which the quark evolves into a jet is called "hadronization" and consists 
of a parton shower, which can be perturbatively calculated, and a fragmentation process, which is a non-perturbative 
process modeled using Monte-Carlo techniques. As will be seen later in this review, the presence of jets in the final 
state gives rise to several experimental challenges. 

Most hadrons are unstable and decay to lighter particles. The probability of decay occurring at a time t is described 
by an exponential, P ~ e~*/ A , where the lifetime, A, is characteristic of the decaying particle. For example, hadrons 
containing a &-quark decay with a lifetime of about 1.5 ps. This is long enough that their flight distance in the lab 
frame is large (of order mm) compared to the experimental resolution - a fact exploited to identify 6-hadrons in jets. 
The ability to identify 6-jets is important in reducing background contributions for many of the tt measurements we 
will discuss. 

From 1989-2009 top quarks could only be produced at the Tevatron, a proton (p) anti-proton (p) collider with 
a center-of-mass energy of ^fs — 1.96 TeV located at Fermilab. For pp collisions at this energy top quarks are 
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predominantly produced in pairs (it) through the strong interaction diagrams shown in Fig. [2j The most recent 
production cross section calculations are accomplished at next-to-leading-order (NLO) and include l ogarithmic cor- 



rections that account for contributions from higher ord er soft gluon radiation ( Ahrens et all . l2010t Cacciari et al 



2008 : Kidonakis and Yost} 2008 : Langenfeld et al . 20091 ). For a top-quark mass of 172 GeV/c 2 the production cross 



section is calculated in iLangenfeld et all l2009f ) to be a(pp -» tt) = 7.57^° pb for y/s = 1.96 TeV. The uncer 



tainty includes contributions from varying the factorization and renormalization scales, and from varying the parton 
distribution functions (PDF). The PDFs describe the momentum distribution of the quarks and gluons that make 
up the protons and anti-protons. They are empirically determined from fits to deep-inelastic-scattering (DIS) cross 
section measurements performed by a variety of experiments. The fits to the DIS data are performed b y a number of 
differ ent collabora tions and are made a vailable as software packages. Common choices are the MRST ([Martin et al 



2007h and CTEQ (jPumplin et al 



20021 ) PDF fits. The cross section above uses the CTEQ6.6 (|Nadolskv et aLl . l2008l ) 
PDF parameterization and changes by approximately =f0.24 pb for a ±1.0 GeV/c 2 change in the top-quark mass. 

In the SM the top-quark branching fraction is completely dominated by t — > W + b with the decays t — > W + q 
(q = s, d) contributing at the < 1% level. Thus for tt events the experimental final state is determined by the decay 
of the W bosons. The leptonic decays (£ — e, r) of the W boson have a branching fractio n of B(W — > Iv) = 10.8% 
each, while the hadronic decays have a total branching fraction of B(W — > qq') — 67.7% ( Nakamura et all l2010h . 
For it events there are then three possible final states. The dilepton final state (dil) corresponds to both W bosons 
decaying leptonically, tt — > i + vW~Vb, and occurs 10% of the time. The lepton+jets final state (Ijt) corresponds to 
one W decaying leptonically and the other hadronically, tt — > Ivqq'bb, and occurs 44% of the time. The all hadronic 
final state {had) corresponds to both W bosons decaying hadronically, tt — > qq'bqq'b, and occurs 46% of the time. 
Figure [3] shows the Feynman diagrams corresponding to these three final states, which, experimentally, are treated 
separately since the contributing sources of background and the dominant detector effects differ among them. Since 
the experiments cannot efficiently reconstruct decays of the tau lepton it should be understood in what follows that 
£ = e, fj. only unless otherwise noted. 

The decay width of the top quark is also well specified in the SM. I ncluding radiative QCD and QE D corrections the 
decay width for a top-quark mass of 172 GeV/c 2 is T t = 1.4 GeV/c 2 (jJezabek and Kfihr] . ll993l[l9M) . The theoretical 
uncertainty is at the level of a few percent and is dominated by uncertainties in a s and in missing higher order QCD 
corrections. The large decay width of the top quark has important experimental consequences. Since T t is large 
relative to the hadronization scale of QCD, Aqcd = 250 MeV, the top quark decays before quark/anti-quark bound 
states are formed. Essentially, the top quark is produced and decays as a free quark - it is unique among the quarks 
in this respect. As a result of this unique feature, many of the experimental techniques used to explore the properties 
of the lighter quarks are not useful in exploring the top quark so that new methodologies had to be developed. 



III. NEW PHYSICS AFFECTING TOP-QUARK PROPERTIES 



It is important to note that the SM has several sh ortcomings an d is widely regarded as a low energy approximation 
of a more complete description of particle physics ( Quigg . 20091) . There are numerous hypothesized theories that 



address these shortcomings and offer this more complete description. The theories differ primarily in the manner by 
which they break the electroweak symmetry and impart mass to the elementary particles. The search to discover and 
describe the physics that accomplishes this is among the most ardently pursued by high energy experimentalists and 
theorists alike. Is it the SM Higgs, or is it some new physics, something Beyond the Standard Model (BSM)? The 
principal aim of top-quark physics is to help answer these questions. With a Yukawa coupling of order unity, the top 
quark strongly couples to the dynamics of electroweak symmetry breaking and is thus expected to be a good place to 
probe for BSM effects. In addition, because of its large mass, top-quark observable can receive large quantum loop 
corrections from possible new particle contributions - significantly larger than for the lighter quarks and leptons. The 
strategy for probing for BSM effects using the top quark employs both direct searches for new particles and indirect 
searches by looking for deviations of experimental observables from the SM expectations. 



New physics contributions can affect top-quark observables directly in two ways (jHill and Parkel 119941) . First non- 
standard top-quark production can appear through intermediate heavy states such as new gauge bosons that decay 
into a tt pair or into a final state mimicking the SM it signature. Secondly the top quark can decay into exotic 
particles such as a charged Higgs boson, t — > H + b, in models with multiple Higgs doublets. In the first case, BSM 
effects distort the SM top-quark production observables, while in the second case BSM effects produce discrepancies 
among the different decay observables. It is thus important to measure the top-quark properties using as many decay 
final states as possible and to compare them. 

There are many examples of new physics models that affect the production or decay observables in tt samples. For 
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example, new heavy states that carry color charge can modify the it production observables. They can be color singlet 
scalar or vector particles or color non-singlet particl es. Examples of models that predict color singlet scalar particles 
are th e supersymmetric extension of the SM, SUSY (jChung et all . 2005t lHaber and Kanei Il98d [Martini . Il997t iNillesl . 
1984), or more generally models with multiple Hig gs fields, MHDM (Grpssmanl.ll994h. Example s of models that predic t 



color singlet vector particle s include te chnicolor (jEichten et al . 19861: Eichten and Lane . 198Clt ). top color drlill . 1991 ) 



or top assisted technicolor dHilll. I1995T). Color non-singlet resonances arise in these BSM theories (|Chivukula et al. 
' 120071: ISimmoni 



19961 : IChoudhurv et al 



1997). Compactified extra dimension models can also affect tt production 



and decay. These models produce Kaluza-Klein (KK) modes that either decay preferentially into tt pairs or into final 
states that mimic the SM s ignature. The fiat extra dimension TeV -1 models produce KK modes that can decay 



ry 

1994 Il999t ICheung and Landsbergl 120021 : iDienes et all Il998l . Il999t iRizzo and WellsL 

tj so that they 



to tt pairs (|Antoniadis et al 

20001) . Universal extra dimension models pair produce a KK mode that ea ch decay as KK 



2, 
(] 



2001; Cheng et 



mimic the SM sign ature but also include additional photons in the final state (jAppelquist et al. 
20021 IrIzzoI. l200ll) . Warped Randall-Sundrum extra dimension models predict KK gravitons that can decay to it 
Davoudiasl et aZ.f . 12001 : Randall and Sundruml . 1999). while bulk Randall-Sundrum models with fermions and gauge 
boso ns in the bulk predict in addition KK top quarks that would be produced in pairs and would decay to Wb as in the 



SM dAgashe et all, 20031: IChang etaH l2000l;ICsaki et all, l2002HDavoudiasl et all l2000l iGrossman and Neubertl . l2000l: 



Hewett , et all bOOa iHuber and Shafil . I2OOII: (Pomaroll l2000l: iRandall and Schwartzl l200lL l2002h. So me BSM theories 



predict a fourth family of (heavy) quarks and leptons (lFrampton et al 



2000; 

rr — 



Hold om et al. 

rr — 



can still be compatible with electroweak precision tests (lEberh a rdt et all 12010; Krib s et al 



2009j). Such extensions 
200 it ). These additional 



generations can appear for instance in grand-unified theories (jLangacken . 119811 ) but some of the models cited above 
also pre dict new particles that look like heayier qu arks. Similarly, Little Higgs models predict new heavy vector-like 
quarks ( Dobrescu and Hill . 19981 : Han et al. , 2003 ). All of these can produce decay final states that mimic the SM it 
signatures. 

New physics contributions can also affect top-quark decays or other prope rties of the i t sample. For example , 
some SUSY and MHDM theories produce top-quark decays into new particles (jChung et all l2005t iGrossma nl. ll994l) . 
Other models modify the SM top-quark branching fractions. For example flavor changing neutral current (FCNC) 
decays, t — > cX°, or t —> uX° where X° is any neutral particle like a photon or a Z boson, can be modified by 
BSM physics. Since these processes are extr emely suppressed i n the SM even small new physics effects can yield 
large changes in FCNC branching fractions ( Aguilar-Saavedral 2004 ). The it sample can also be probed for the 
presence of more general anomalous couplings ( Kane et al\ , 1992 ) affecting the kinematic properties of top-quark 
decays and appearing as discrepancies with the SM predictions. These might affect the it forward- backward charge 



asymmetry, the tt spin correlations , or the polarization fra ctions of W bosons from top-quark decays ( Cao et al. , 2010t 



Cheung et all 1200 



Diouadi et all 120091 : IShu et aLl . l2010h . The measured intrinsic properties of the top quark might 
also be affe cted. For example some BSM theories predict particles with electric charge —4/3 or +5/3 that decay into 
top quarks (jContino and Servantl . 120081 ) . 

Over the last decade the Tevatron experiments have undertaken a broad top physics program aimed at exploring 
these possibilities both by comparing measured top-quark properties to SM predictions, as discussed in Sec. IIX( and 
by directly searching for new particles in the top-quark event samples, as discussed in Sec. Ixl 



IV. EXPERIMENTAL APPARATUS 

The pp collision data at the Tevatron are collected and analyzed by two collaborations, CDF and D0. For both 
the CDF and D0 experiments the detector apparatus is designed to be a general purpose collider detector with the 
capability of reconstructing charged particle trajectories and jets, identifying electrons, muons, photons, and 6-quark 
jets, and measuring the energy and momentum of as many of the final state particles as possible. Functionally the 
two detectors behave very similarly, although they differ significantly in their technical details. They both have a 
cylindrical geometry centered on the pp beam line and employ a coordinate system with the z axis parallel to the beam 
line and pointing in the direction of the proton beam. The angle 9 is the polar angle relative to the z axis and the angle 
(j) is the azimuthal angle in the plane transverse to the z axis. The energy, E, and momenta, p, of electrons, muons, 
photons, and jets are estimated using information the calorimeters and spectrometers. It is convenient to define the 
transverse energy, Et — Esin9, and the transverse momentum, px — \p\sin9. Since the pp beams are traveling 
along the z axis the initial state is known to satisfy (Et,pt) — (0,0) while the (E z ,p z ) components of the initial 
state are unknown cvent-by-event - and in general are non-zero - due to the ambiguities introduced by the parton 
distribution functions. Because the collision events are boosted along the z-axis hadron collider experiments often use 
pseudo-rapidity, 77 = — In tan |, in place of cos 9 when defining angular regions since differences in pseudo-rapidity, 
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A?7, are invariant under these boosts. 



T he CDF and D0 dete ctors are shown in Fig. 0] and are described in detail in references (jBlaire et all 119961 ) 



and ( Abazov et all l2006aj) and only a generic description is provided here. The detectors are divided into a central 
region covering (approximately) |ry| < 1.3 and a forward region covering (approximately) 1.3 < |r/| < 3.5. In the 
central region charged particle trajectories are reconstructed as tracks using a set of custom built components. At 
the inner most radii, immediately surrounding the pp beam line, are precision vertex systems comprised of several 
concentric layers of silicon microstrip detectors with a per layer charged particle position resolution of around 10 /im 
in the transverse plane. These detectors are important in reconstructing the primary pp interaction vertex as well as 
identifying displaced vertices from the decay of long lived particles (e.g. &-hadrons). The silicon systems are surrounded 
by large radii, multi-layered, low mass tracking chambers which achieve position resolutions of around 150 /im per hit. 
The silicon and tracking systems are inside superconducting solenoids, which provide magnetic fields parallel to the 
beam line and induce curvature in the trajectories of charged particles from which their momentum can be deduced. 
The transverse momentum of charged particles is typically measured with a resolution of a PT /pT < 0.01 px- Outside 
the solenoid are the calorimeters with an electromagnetic compartment in front of a hadronic compartment. The 
electromagnetic (EM) compartment consists of roughly 20 radiation lengths or more of material and fully contains 
e^- and 7 showers. The showers induced by pions, kaons, and other long lived hadrons are initiated in the EM 
calorimeter but continue through and are contained in the hadronic compartment. The energies of electrons and 
photons are typically measured with a resolution of approximately cr E /E = 14%/y/E while the energies of jets are 
typically measured with a resolution of approximately <je/E = 10 — 15%/y/E. The calorimeters are segmented into 
towers with a typical granularity of rj x <p = 0.1 x 0.1. Beyond the hadronic calorimeters are the muon chambers, 
small scintillator plus tracking systems which serve to identify muons, the only charged particles able to traverse 
the calorimeters without being absorbed. In the forward region there is some limited tracking capability from the 
silicon system, but with a much worse px resolution than that achieved in the central region. There are also EM and 
hadronic calorimeters together with muon chambers. In the far forward region, \rj\ > 3.5, is a luminosity monitor used 
to measure the pp collision rate in the detector during data taking. In total each experiment has about one million 
read-out channels. 

The data acquisition system digitizes the read-out channels and records them to permanent storage (tape) for 
later analysis. The Tevatron delivers a collision to each experiment every 396 ns, which corresponds to about 2.5M 
collisions per second. This is far too many collisions to record them all and each experiment employs a three level 
trigger system to identify the most interesting events to keep for later analysis. The level one trigger is hardware 
based and employs simple thresholds and a coarse granularity to identify events with high energy depositions in the 
calorimeter or segments in the tracking system consistent with having originated from a high momentum particle. The 
level one system reduces the initial event rate by about two orders of magnitude. The level two trigger system employs 
some custom hardware and processing to match objects identified at level one and to achieve improved resolution 
with finer granularity. For example, at level two a high momentum track may be matched to a high energy deposition 
in the EM calorimeter or to a track segment in the muon chamber, consistent with what is expected from an electron 
or muon, respectively. The level two system reduces the event rate by another one order of magnitude. The level 
three system employs a computing farm and custom software to reconstruct the full event and confirm the level two 
decision with a resolution and granularity that nearly matches that achieved after final calibrations and offline event 
reconstruction. The level three system achieves an additional order of magnitude reduction in the event rate so that 
each experiment writes the most interesting events to tape at a rate of about 100 Hz. 



V. EVENT RECONSTRUCTION 

The events which are written to tape are later processed on large scale computing farms using custom software 
to perform event reconstruction. This offline processing uses final calibration constants and alignments and more 
sophisticated algorithms relative to what is used in the level three trigger in order to achieve improved efficiencies and 
resolutions. The reconstruction software combines information from the tracking, calorimeter, and muon chamber 
systems to provide particle identification. Charged particles reconstructed in the tracking systems are extrapolated 
to the calorimeters and muon chambers. Tracks matched to energy depositions in the EM calorimeter or to muon 
chamber hits are taken to be electron or muon candidates, respectively. An EM energy deposition without a track 
match is taken to be a photon candidate. A contiguous set of calorimeter towers with energy depositions above a 
predetermined noise threshold and matched to a collimated set of tracks in the tracking chamber is taken to be a jet 
originating from the hadronization of an initial quark or gluon state. A jet which contains tracks that form a secondary 
vertex significantly displaced from the primary interaction vertex is taken to be a candidate 6-jet. The presence of 
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a neutrino in the final state is inferred from a momentum imbalance in the transverse plane since neutrinos traverse 
the entire detector volume without interacting. In the events of interest the final state particles are all relativistic, 
so the momentum imbalance is calculated using the calorimeter energies and is called "missing Et" or MET and is 
defined as BJt = — J2i E T i, where the index i runs over all calorimeter towers, E T is the transverse energy measured 
in the ith tower and i is the unit vector pointing from the primary event vertex to the center of the ith tower in the 
transverse plane. The resolution of the E^t calculation can be improved by correcting for muons and jets in the event. 

The it all hadronic final state (it — > qq'bqq'b) is characterized by a large jet multiplicity, and the presence of 2 
6-quark jets. Analyses in this final state exploit the jet reconstruction and 6-jet identification capabilities of the 
detector to identify events of interest. The it lepton+jets (it — > ivqq'bb) final state contains one high energy lepton, 
large E<t> high jet multiplicity, and 2 6-quark jets. The it dilepton (it — > t + vbl~Vb) final state has two high energy 
isolated lcptons and large E^t together with 2 6-quark jets. Analyses using the dilepton or lepton+jets final state 
exploit the full range of detector capabilities using electron, muon, and 6-jet identification, jet reconstruction, and 
missing energy estimates to identify candidate it events. Each analysis begins by choosing the subset of recorded 
data satisfying a particular set of trigger criteria. This eliminates the need to run the full analysis software over those 
events whose topologies and kinematics significantly differ from those expected of it events. The trigger criteria are 
designed to exploit a given characteristic of at least one of the it final states and to identify those events with high 
efficiency. For each event satisfying the trigger criteria the primary interaction vertex is identified, leptons and jets 
are reconstructed, an estimate of E^t is made, and 6-jets are identified using a dedicated algorithm. The event sample 
at this stage is dominated by background events originating from non-it processes. Additional final state dependent 
selection criteria are used to suppress these background contributions and improve the purity of the sample. 



A. Trigger Criteria 

The had data set is collected using trigger criteria requiring at least four jets each with raw (uncorrected) energy of 
Et > 15 GeV in events with large total transverse energy as determined by summing over all calorimeters towers (e.g. 
^2iE T > 175 GeV). For each 1 pb _1 of collision data these triggers collect about 8-10 thousand events, depending 
on the details of the trigger selection, dominated by QCD production of multi-jet final states events (e.g. pp — > qqqq 
or pp-> qqgg). 

The dil and Ijt data sets are collected using several different sets of trigger criteria. The two most important 
sets of trigger criteria focus on identifying the electron or muon from the IV-boson decays^. These leptons are 
typically high energy and isolated. The first set of trigger criteria is an inclusive high-i/r electron data set requiring 
a candidate electron with Et > 20 GeV. The electron is identified at the trigger level as a high energy cluster in the 
EM calorimeter matched to a high momentum track in the tracking chamber. Additional requirements are made on 
the ratio of energy deposition in the EM calorimeter to the energy deposition in the hadronic calorimeter, and on the 
shape of the EM shower in the calorimeter. The electron candidate must be isolated, with few surrounding tracks 
and/or only small additional energy deposits in neighboring calorimeter towers. The second set of trigger criteria is 
an inclusive high-p-r muon data set requiring a candidate muon with pt > 20 GeV/c. The muon is identified at the 
trigger level as a track segment reconstructed in the muon chambers matched to a high momentum track reconstructed 
in the tracking chamber. The muon candidate is additionally required to be isolated. Because these first two sets of 
trigger criteria make tracking requirements their acceptance is limited to the central detector region. The number of 
it events collected can be significantly increased by including triggers which collect events in the forward detectors. 
In order to keep the purity of the trigger sample reasonable these forward triggers are less inclusive. They usually 
require some high energy jets (e.g. Ej? > 10 GeV) in addition to a high energy EM calorimeter cluster in the forward 
region (e.g. E^ M > 20 GeV) or large Ufa (e.g. S^t > 35 GeV). For each 1 pb _1 of collision data these triggers 
collect approximately 1-2 thousand pp — > W + jets and pp — > Z/j + jets events, depending on the details of the trigger 
selection. 



Because of the additional experimental challenges they introduce, it events with one of the W bosons decaying to a tau lepton are not 
often used. Since the tau lepton quickly decays to a multi-particle final state the W —y tu decay is difficult to trigger on in an inclusive 
manner. Those leptonic tau decays, r — > evu and r — ► fj,uu, satisfying the inclusive muon and electron trigger criteria are included in 
the data set described here. 
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B. Primary Interaction Vertex 

The counter-rotating proton and anti-proton beams at the Tevatron are brought into collision near the centers of 
the CDF and D0 detectors. The exact location of the collision varies for each beam crossing. The shape of the 
luminous region is determined using tracks in di-jet events (pp — > qq), which are predominantly promptly produced 
and thus originate from the collision vertex. The tracks are constrained to a common origin and the distribution of 
the (x, y, z) position of the resulting vertex is determined using millions of eventt0. The distribution of the collision 
positions along the z axis can be approximated as a Gaussian with a mean near z = and an RMS of approximately 
30 cm. In order to achieve the best jet energy resolutions and highest &-jet identification efficiency it is necessary to 
reconstruct the location of the interaction vertex for each event. For dil and Ijt events the z coordinate of the highest 
energy lepton track can be used together with the known shape of the luminous region, as determined from the di-jet 
data as described above, to estimate the location of the interaction vertex with a resolution of a few hundred microns 
in z and about 35 \xm in x and y. These resolutions are improved by a factor of three to five by using the tracks in 
the event itself. In order to avoid bias from long lived 6-hadron decays in &-jets, only tracks whose distance-of-closest- 
approach (also known as the impact parameter, do) to the beamline is comparable to the do resolution are used. These 
tracks are then constrained to originate from a common point, whose coordinates are taken as the position of the it 
interaction vertex. An iterative procedure is employed by eliminating tracks with large do relative to the fitted vertex 
to form a new vertex. The process is repeated until no large impact parameter tracks are included in the vertex fit. 
For had events there is no lepton to seed the identification of the it vertex. Again, tracks in the event are used and 
are constrained to originate from a common origin and an iterative procedure is used to remove the bias from tracks 
originating from 6-hadron decay. Depending on event topology, resolutions of 5 — 20 /Ltm in the transverse plane and 
of 100 /im along z are achieved on an event- by- event basis. 

Since the pp inelastic cross section is very large there is often more than one hard scatter interaction per beam-beam 
crossing. The mean number of expected hard scatters is linearly dependent on the instantaneous luminosity of the 
colliding beams. For the Tevatron parameters the mean number of expected hard scatters per beam-beam crossing 
is about 3.5 at an instantaneous luminosity of A ns t = 1.0 x 10 32 cm _2 s _1 . The actual number of hard scatters 
per event is Poisson distributed about the mean. These additional interactions are typically separated in z. Since 
the RMS of the luminous region along the z axis is large relative to the resolution with which a single interaction 
vertex is reconstructed, the experiments can efficiently reconstruct and readily resolve multiple interaction vertices 
in a single event. This is important since the typical Tevatron stores have an initial instantaneous luminosity that 
exceed 3 x 10 32 cm _2 s _1 and an average instantaneous luminosity over the duration of a store, which is typically about 
15 hours, of about 1.5 x 10 32 cm _2 s _1 . For events with multiple reconstructed interaction vertices it is necessary to 
identify which of the vertices corresponds to the it interaction vertex. For the leptonic final states the vertex closest 
to the high energy lepton is used. For the had final state, the scalar sum pt of all tracks associated with each vertex 
is calculated and the vertex with the largest ^ \pt\ ls use d- The it interaction vertex is correctly identified nearly 
100% of the time. 

C. Electrons and Muons 

The reconstruction of high energy electrons and muons from VT^-boson decays is important in identifying the signal 
it events. The leptons from Z e + e~ and Z — > fi + pT decays are kinematically very similar to the W decay leptons 
and are thus used as data control samples to measure efficiencies and backgrounds for the lepton reconstruction 
and selection criteria employed in the it analyses. The reconstruction of lower energy electrons and muons from 
semi-leptonic 6-hadron decays is sometimes used to identify 6-jets. 

Electron reconstruction begins by identifying energy depositions in the EM calorimeter above a specified threshold, 
typically around 2 GeV. In the central detector region, tracks from the tracking chamber with momentum above a 
few GeV/c are extrapolated from the tracker outer radius to the front face of the EM calorimeter. If the extrapolated 
track position is consistent with the energy-weighted centroid of the calorimeter cluster within uncertainties, including 
contributions from multiple scattering in the material between the tracker and the calorimeter, the resulting EM 
calorimeter and track match are taken to be a candidate electron. To reduce contributions from hadrons mimicking 
electrons, the shape of the shower in the calorimeter, and the ratio of track pt to calorimeter Et are required to be 
consistent with that expected of electrons as determined from Z — > e + e~ events and test beam data. The energy 



2 For a vector of magnitude v we define (v x , v y , v z ) = t>(sin cos 0, sin0sin0, cos0). 
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deposited in the hadronic calorimeter tower immediately behind the candidate electron is required to be small relative 
to the EM energy deposit. Electrons consistent with having originated from 7 — > e + e~ conversions in the detector 
material are identified with dedicated algorithms and removed from further consideration. In the forward detector 
region where there is very limited coverage from the full tracking chamber, electrons are identified using stricter criteria 
on the EM cluster energy, shower shape, and hadronic to EM energy deposition ratios. The purity of the resulting 
forward electron sample can be significantly improved by searching for hits in the tracking system consistent with 
expectations derived from the Et and position of the EM cluster and the location of the event primary interaction 
vertex. 

Muon reconstruction begins by identifying track segments in the muon chambers. Tracks from the tracking chamber 
with momentum above several GeV/c are extrapolated from the tracker outer radius to the front face of the muon 
chamber. If the extrapolated track position is consistent with the position of the muon chamber track segment 
within uncertainties, including contributions from multiple scattering in the material between the tracker and the 
muon chamber, the resulting segment and track match are taken to be a candidate muon. Through going cosmic 
ray muons are identified using dedicated algorithms and removed. Hadrons whose shower are not fully contained in 
the calorimeter and other intervening material are said to "punch-through" to the muon chambers and can mimic a 
muon signature. To reduce contributions from punch-through hadrons, stricter segment-track matching and track pr 
criteria together with track quality criteria can be employed. The track quality criteria usually include requirements 
on the number hits used to form the track, the impact parameter, and the track fit \ 2 ■ These track quality criteria also 
help reduce contributions from hadron decays-in-flight (e.g. K — > n^v). The energy depositions in the calorimeter 
towers traversed by the muon can be required to be consistent with that expected from a muon as determined using 
Z — > fi + ji" and J/tp — » data control samples. 

The leptons from W- and Z-boson decays are typically well separated from other objects in the event. To further 
suppress hadronic backgrounds isolation criteria are employed. Calorimeter isolation requires that the energy deposited 
in an annulus around the calorimeter towers associated with the lepton candidate be small. Tracking isolation requires 
that the scalar sum pt of all tracks in a cone around the lepton candidate be small. The requirements are sometimes 
absolute (e.g. -E(annulus) < 3 GeV) and sometimes relative (e.g. i?(annulus)/-E(electron) < 0.1). 

Samples of Z —> l + t~ events collected with the inclusive electron and muon triggers described above, are used to 
measure the trigger, reconstruction, and identification efficiencies for electrons and muons, respectively. The "tag- 
and-probe" method is employed where the lepton leg which fired the inclusive trigger is the "tag" leg, while the 
second leg is unbiased and used as a "probe" to measure the efficiency of interest. The efficiencies are measured 
with total uncertainties of < 1% relative and are characterized as a function of pr, T), instantaneous luminosity, 
an d number of reconstruct ed interaction vertices. A more detailed description of these measurements can be found 



in (jAbulencia et all l2007bf ) 



The absolute momentum scale of the spectrometer is determined using J/ip, T, T', and Z-boson decays to 
by comparing the reconstructed invariant mass of the muon pairs to the known masses of each resonance. By 
using a variety of resonances any non-linearities in the momentum scale can be identified and corrected for. These 
studies result in a more precise understanding of the magnetic field uniformity, the distribution of material in the 
detector volume, and any small residual mis-alignments in the tracking systems. The absolute energy scale of the 
EM calorimeter is determined using ir° —> 77, J/ip —> e + e~ and Z —> e + e~ decays and using the E/p distribution of 
electrons from W — ¥ ev decays. Again, by using a variety of decay processes non-linearities in the energy scale can 
be identified and corrected for. In both cases the absolute scales are determined to a precision well below 1% and are 
monitored continuously and kept constant at the 0.2 — 0.3% level. The observed RMS widths of the reconstructed mass 
distributions are dominated by the detector resolutions. Studies of the di-muon and di-electron resonance widths help 
determine the momentum and energy resolutions, respectively. The absolute scales and resolutions are characterized 
as a function of energy /momentum, 77, instantaneou s luminosity, and number of reconstructed in teraction vertices. 



These methodologies are described in more detail in (jAaltonen et all l2008ct lAbazov et all l2009al) . 



D. Jets 

Jets are reconstructed as a contiguous set of calorimeter towers each above a predetermined noise threshold. A 
variety of algorithms are available to determine which towers to group together. The algorithms differ in the criteria 
they use to combine towers to form jets and in the manner in which they identify and merge or separate jets which 
are near each other in 77 — 4> space. The behavior of the algorithms, particularly with regard to their handling of 
nearby jets, affects the degree to which they are sensitive to low energy "soft" QCD radiation. Algorithms which 
are insensitive to soft QCD radiation are said to be "infrared safe" and can be compared to theoretical predictions 
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in a very clean manner. There are in general two classes of jet algorithms. The kx algorithms combine calorimeter 
towers to form jets based on their relative transverse momentum as well as their rj — (j) separation. They are rigorously 
infrared safe. The cone algorithms combine calorimeter towers to form jets based only on their 77 — cf> separation. 
The cone algorithms begin by identifying a set of seed towers with energy above some seed threshold, typically about 
1 GeV. The Midpoint algorithm employed by CDF and D0 is such an algorithm with additional steps introduced in 
order to minimize the sensitivity to soft QCD radiation. While rigorously the Midpoint algorithm is not completely 
infrared safe, in practice the residual sensitivity is experimentally inconsequential. This is due to the fact that the 
energy of jets in tt events is large (approximately 50 GeV on average) compared to the seed threshold so that small 
fluctuations in tower-to-tower energy depositions from soft QCD radiation do not affect the results of the algorithm. 
The cone algorithms offer an experimental advantage because the physical extent of each jet is well defined to have a 

radius of R = y (Ar/) 2 + while the physical extent of jets formed using a for algorithm vary from jet to jet. In 



practice the two algorithms yield very similar results ( Aaltonen et ai , 2008f). In the cone algorithms the cone radius 



employed can be varied. Since the jets in tt events are typically collimated a relatively narrow cone size is employed, 
R = 0.4 for CDF and 0.5 for D0. 

The jet reconstruction algorithm uses raw calorimeter energies. The total energy of the jet is estimated as the 
sum of energies over all towers included in the jet. The direction of the jet is determined using an energy weighted 
average over all contributing calorimeter towers. The energy and direction of the jet are then used as an estimate of 
the four momenta of the parent parton (ie quark or gluon). The jet energy resolution is improved by employing a set 
of corrections which account for several effects. The effects fall into two broad categories, instrumental effects, and 
physics effects. The instrumental effects accounted for include corrections for non-uniform and non-linear response of 
each calorimeter tower as well as contributions from multiple interactions per event. The physics effects accounted for 
include corrections for contributions from the underlying event, from the hadronization process, a nd from out-of-cone 



show ering. The cal i bration of the CDF and D0 jet energy corrections are described in detail in ([Bhatti et all 120061 ) 



and (jAbazov et all l2008af) and are a function of jet transverse energy and pseudo-rapidity. The corrections vary 



from about 60% for jets with transverse energy around 20 GeV to about 20% for jets with transverse energy around 
100 GeV. The corrections can be checked using pp —> 7+ jet and pp —> Z +jet events. Given the initial constraint that 
(Et,Pt) — (0,0), it is expected that the jet transverse energy after all corrections should be equal to the transverse 
energy of the well measured photon or Z boson (when Z/-f — >• e + e~ or /i + /i~) in the event. The uncertainties on the 
jet energy scale derive from many sources, some of which are experimental (e.g. the statistics of the control sample 
or an inadequate understanding of the detector response) and some of which are modeling (e.g. modeling of the 
underlying event, hadronization, and out-of-cone showering processes) related. Like the corrections, the uncertainties 
are jet Et and r\ dependent and range from about 10% at low energies (< 40 GeV) to a few 1% at high energies 
(> 75 GeV). The dominant source of uncertainty is modeling related at the lower energies and experiment related at 
higher energies. 

The full set of jet energy corrections are designed to give an accurate estimate of the energy and momentum of the 
parent parton that had initiated the jet. However, the full set of corrections are not always employed. Frequently 
jets are only corrected to the particle level, which effectively corrects for all the instrumental effects as well as contri- 
butions from the underlying event and multiple interactions. Jets corrected to this level are effectively independent 
of apparatus and environment and can be compared across experiments. For the selection criteria described below, 
the jets are corrected to the particle level. Later, when using the selected tt candidate events in analyses the jets are 
corrected to the parton level whenever a kinematic fit is employed. 



E. Missing Transverse Energy 

We exploit the fact that the initial state satisfies (Et,Pt) = (0,0) and use an observed imbalance of momentum 
in the transverse plane to infer the presence of high energy neutrinos in an event. As mentioned above, because the 
final state particles are all relativistic at the relevant collider energies, we use the energy measured in the calorimeters 
as an estimate of the total momenta of particles produced (including neutral particles). For this reason the quantity 
examined is called missing transverse energy, S^t- The initial estimate of the missing transverse energy begins with 
using only calorimeter information: B}t = — 5Zj E T i, where the index i runs over all calorimeter towers, E T is 
the transverse energy measured in the ith tower and i is the unit vector pointing from the primary event vertex 
to the center of the ith tower in the transverse plane. Corrections are made if there are identified muons in the 
event since they are minimum ionizing and the calorimeter energy deposition is a poor estimate of their momenta: 

BJt = ~(^2iE T i —E^i^+p T i^), where E T , i^, and p T are the calorimeter transverse energy, azimuthal angle 
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unit vector, and tracker transverse momentum associated with the muon candidate. The Bjr resolution can be 
further improved when jets have been reconstructed by taking advantage of the jet energy corrections described 
above by removing from the initial sum any calorimeter tower included in a jet and substituting the jet transverse 
energy (corrected to the particle level) at the jet <j>. For tt analyses the final expression is often rewritten: S^j, = 

— wZiPt^ + Eip Cm ij + J2uc -E"r c *uc^ j where £ is the sum over identified leptons (electrons or muons) in the event, 
j is the sum over reconstructed jets, and uc is the sum over unclustered calorimeter towers (ie. those not included 
in a jet). For electrons (muons), p T is estimated from the EM calorimeter (spectrometer). For jets, E^ C01 is the jet 
transverse energy corrected to the particle level. For the unclustered calorimeter towers the raw calorimeter energy 
is used. This expression is equivalent to the original expression with corrections for the muons and jets as described. 

At the trigger level only the raw S^t estimate is available, without any of the corrections. At the event selection 
and later in the analyses employing kinematic fits, the fully correctly B}r is used. 



F. 6-quark Jets 

Most jets are initiated by a parent light quark (uds) or gluon (g). Thus the identification of at least one of the 
6-jets in tt events significantly reduces contributions from background processes. The 6-jet identification algorithms 
employed by the Tevatron experiments generally exploit two main features of 6-hadron decay - their long lifetime 
and their semi-leptonic decays. While c-hadrons also have some of these same features, on average their lifetimes are 
shorter and they contain fewer electrons and muons than 6-hadrons so that the algorithms employed typically achieve 
6-identification efficiencies a factor of 2-4 larger than the corresponding c-identification efficiencies, which is sufficient 
to reduce charm quark related backgrounds to acceptable levels for the tt analyses described here. 

The branching fraction for a 6-hadron into e or fi plus anything is about 20% each (including sequential 6 — > c —¥ eX 
or (J.X), which is about twice as large as that for c-hadrons, while udsg-jets have a negligible e/fj, content. The 
energy of these leptons from 6-decay is a few GeV on average, so much "softer" than those from W- or Z-boson 
decays. Consequently these algorithms are often called "Soft Lepton Tagging" algorithms. Additionally, since they 
are embedded in jets, these leptons are less isolated. The identification of these electrons and muons uses the same 
reconstruction algorithms as those discussed in Sec. IV. CI To further suppress backgrounds from hadrons mimicking 
electron or muon signatures these algorithms demand the leptons to have a transverse momentum greater than about 
3 GeV/c and make stricter matching requirements between the track and EM calorimeter cluster (electrons) or track 
segment in the muon chamber (muons) . In addition the EM shower shape is exploited to suppress backgrounds for the 
electron case. The full set of identification criteria is usually combined in a single variable using a relative likelihood 
function, a global \ 2 ° r other similar multivariate techniques. Data control samples of J/ip — > fi + /j,~, 7 — > e + e~ , 
and Z — > /J, + /i~ , e + e~ are used to measure the identification efficiencies. Typical efficiencies achieved are 80-90% for 
muons in 6-hadron jets above the px threshold and fiducial to the muon chambers. The efficiencies for electrons in 
6-hadron jets is about a factor of two lower because the EM tower often includes contributions from nearby hadrons 
in the jet thus degrading the shower shape criteria. The degree to which this is a problem will partly depend on the 
segmentation of the EM calorimeter. For both algorithms the rate at which hadrons are mis-identified is about 0.5% 
per track and varies by about ±20% (relative) as a function of track p?, track 77, or hadron species (ie. ir^ 1 vs K ± 
vs p/p). These fake rates are determined using samples of A — > pir~ , D + * — > D q tt + — » K~ir + ir + , and K ( J — > ir + ir~ 
decays. Including the semi-leptonic branching fractions and the track multiplicities in udsg-jets these algorithms in 
combination achieve an efficiency of ~ 25% per fiducial 6-jet in tt events and misidentify a udsg-\et at a rate o f 



0. 02 — 0.03 each. These algorithms and their performance are described in more detail in (jAcosta et all l2005al ) 



and (jAaltonen et all l2010ff ) 



The most important 6-jet identification algorithms exploit lifetime information by using tracks originating from 
the 6-hadron decay and reconstructing a secondary vertex significantly displaced from the primary interaction vertex. 
These algorithms exploit almost all 6-decay topologies and thus make use of nearly the full branching fraction. They 
work on a per jet basis and begin by associating to each jet a set of tracks using the r\ — 4> separation between the track 
and the jet axis. Only tracks within a given radius (e.g. R < 0.4), above a given pt threshold (e.g. pr > 1 GeV/c), 
and whose z coordinate at the distance of closest approach to the beam line is consistent with having come from the 
primary interaction vertex of the event are used to form the candidate secondary vertices. These tracks are vetoed if 
they form a good vertex and yield an invariant mass consistent with a K® or A when paired with any other track in 
the same jet. They are also vetoed if they are identified as originating from a photon conversion, 7 — » e + e _ , in the 
detector material using dedicated algorithms. At hadron colliders "build-up" algorithms are preferred to "tear-down" 
algorithms because non-6-decay tracks often out number tracks from 6-decay owing to contributions from the under 
lying event and hadronization process. The build-up algorithms begin by taking the surviving tracks pair wise and 
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making a seed vertex, often times beginning with those tracks having the largest impact parameters relative to the 
primary interaction vertex of the event. Additional tracks are then added to these seed vertices. The px and impact 
parameter requirements for these additional tracks may be relaxed relative to the requirements made on the tracks 
forming the seed vertices. Only vertices with good fit x 2 are kept and tracks making large contributions to the \ 2 
are removed. Surviving vertices are then subject to additional selection criteria optimized to suppress the rate of fake 
tags from udsg-jets while maintaining high efficiency for &-jets. In particular the decay length, which is the magnitude 
of the vector connecting the primary interaction vertex of the event to the candidate secondary vertex is required to 



fall within the beam pipe (e.g. 



L 



< 1.5 cm at CDF and 



L 



< 2.6 cm at D0) in order to eliminate interactions 



in the detector material, the angle between L and the jet momentum must be less than 7r/2, and the decay length 

significance must exceed some minimum (e.g. L /ctl > 7). There can be multiple seed vertices per jet, which are 

ordered by vertex quality. Each seed vertex is then tried in turn until a good vertex is found. The algorithms may 
exploit additional properties of the candidate secondary vertex, like the number of tracks included in the vertex or 
the invariant mass of the tracks included in the vertex, to further discriminate 6-decay vertices from backgrounds. 
The efficiency of these algorithms is determined from the data using control samples of 6-jets collected with triggers 
which identify the soft electrons or muons from semi-leptonic 6-decays. These samples thus yield the efficiency for 
semi-leptonic 6-decays and tend to have a fairly soft spectrum. Monte Carlo simulation is used to extrapolate 
the efficiencies to generic 6-decays and higher jet energies. The Monte Carlo and detector simulation are tuned to 
accurately reproduce the behavior of the algorithm as measured in the data control sample. Systematic uncertainties 
account for discrepancies between the data and the Monte Carlo modeling and for the uncertainties associated with the 
extrapolation to higher jet energies. The rate at which udsg-jets are mis-identified by these algorithms is measured 
in the data using a sample of generic QCD jets collected using only calorimetric trigger information. To get an 
estimate of the mis-identification rate the experiments exploit the fact that the observed decay length distribution for 
udsg-jets is nominally symmetric about zero with a width dominated by resolution effects. The negative decay lengths 
correspond to those jets for which the angle between L and the jet momentum is larger than 7r/2 and are dominated 
by contributions from udsg-jets. The rate of negative decay lengths is thus used to parameterize the mis-identification 
rate for udsg-jets reconstructed with positive decay lengths. The mis-identification rate is parameterized as a function 
of several variables, such as the number of tracks in the jet, the jet energy, the jet pseudo-rapidity, the number of 
additional interactions in the events, etc. Small corrections are made to account for the 6-jet and c-jet contributions 
to the negative decay length sample and for contributions to the positive decay length rate from residual photon 
conversions in the material. These algorithms are tuned to several operating points. The most popular operating 
point is labeled "Tight" and achieves an efficiency of ~ 50% per fiducial 6-jet in ti events and mis-identify udsg-jets 
at a rate of < 0.01 each. The efficiency has a total uncertainty of about 7% (relative) dominated by systematic 
contributions from uncertainties in the sample composition of the data control sample, the Monte Carlo modeling, 
and the extrapolation in E^ 1 . The mis- identification rates have a total uncertainty of 10 — 20% (relative) dominated 
by systematic contributions from uncertainties in the sample composition, corrections accounting for trigger biases, 
and corrections for the 6-jet and 7 — > e + e ~ contributions. The se algorithms and their performance are described in 



more detail in (|Acosta et all |2005bJ) and (|Abazov et all l2010bf ) 



Other algorithms exploiting the long lifetime of 6-decays are sometimes employed such at "jet probability" algo- 
rithms, or simple "counting" (displaced tracks) algorithms. In general these offer additional efficiency at the cost of 
higher udsg-jet mis-identification rates. 

Often multiple algorithms are employed and then combined using various multivariate techniques to yield another 
10-20% (relative) increas e in the efficiency for identifying &-jets for the same udsg misidentification rate. An example 
is described in detail in ( Abazov et "q*l l2010bh . 



G. Kinematic Fits 



In many analyses kinematic fits are employed to fully reconstruct the ft kinematics and improve the resolution on 
the reconstructed momenta of the final state partons from the top-quark decays. The fits constrain the kinematics 
of each event assuming the ti hypothesis and using the decay relevant for the final state in question. In the had 
channel the six highest Et jets are assumed to correspond to the quarks from the top-quark and VF-boson decays. 
In the Ijt channel the four highest Et jets are assumed to correspond to the quarks from the top-quark decays and 
the hadronically decaying W boson, while the charged lepton is assumed to come from the leptonically decaying W 
boson. In the dil channel the two highest Et jets are assumed to correspond to the 6-quarks jets from the top-quark 
decays and the two (oppositely) charged leptons are assumed to come from the two leptonically decaying W bosons. 
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In the Ijt and dil channels the measured S}t is used to constrain the x and y components of the neutrino momenta. 
The measured four-vectors for the jets are corrected using the jet energy corrections described in Sec. IV. Dl and the 
four-vectors of the leptons incorporate the calibrations described in Sec. IV.CI The calibration studies also yield 
estimates of the corresponding resolutions. The kinematic fit proceeds via a x 2 minimization. As an example, the \ 2 
function employed in the Ijt channel, it — > ivqq'bb, is 



Similar expressions are used for the had and dil channels except that the and triMv terms are replaced by rrijj 
and rribjj terms or vice versa, respectively. The first term constrains the pt of the lepton and jets to their measured 
values within their assigned resolutions, o^. The second term makes the same constraints on the x and y components 
of the unclustered energy in the event. The S^t is thus derived using these first two terms as discussed in Sec. IV. El 
The last four terms constrain the reconstructed invariant mass of some combination of final state partons to the pole 
mas^H of the object from which they are assumed to have originated. For example, rrijj is the invariant mass of 
the two jets assumed to have originated from the hadronically decaying W boson and so it is constrained to equal 
raw within some tolerance, taken to be the natural width of the W, Tw- Similarly for the decay products of the 
leptonically decaying W boson and the decay products of each top quark. In the fit, the p^, fit and p" c ' fit terms arc 
varied to minimize the x 2 ■ For most analyses the mw and to* values are fixed to the world average measured values. 
In analyses measuring the top-quark mass, discussed in Sec. IIX.C.1I the nit term is left a free parameter, m\ cco , and 
corresponds to the reconstructed top-quark mass most consistent with the kinematics of that particular event. 

In order to perform the minimization it is necessary to assign each jet as having originated from a particular 
final state quark. Since the flavors of the parent quarks are unknown, all possible combinations are tried and the 
combination yielding the lowest fit x 2 is used. This procedure yields the correct jet-parton assignment 65-85% of the 
time depending on the details of the selection criteria and decay channel. This ambiguity in the jet-parton assignments 
gives rise to a combinatoric background from it events for which the incorrect combination is used. The combinatoric 
background significantly dilutes the resolution with which the it kinematics can be reconstructed. The number of 
possible jet-parton assignments for an event with Nj jets is iV com b = Nj\. Since the x 2 does not distinguish between 
a swap in the quark assignments for the jets assumed to have come from the hadronic VF-boson decay, the effective 
number of jet-parton combinations is reduced by a factor of two for each W — ¥ jj decay in the final state. The number 
of combinations considered can be further reduced by rejecting combinations for which a jet identified as a fe-jet is 
assigned as one of the light-quark jets from the hadronic VF-boson decay. For example, in a Ijt event reconstructed 
with 4 jets in the final state, there are ^ = 12 jet-parton assignments considered; this is reduced to six different 
assignments if one of the jets is identified as a 6-jet and only two possible assignments if two jets are identified as 
6-jets. The combinatoric background is less of a problem in the dil channel and is severe in the had channel. 

The presence of neutrinos in the final state introduces a further ambiguity. In the Ijt channel the x and y components 
are constrained using the estimated S^t, but the z component is unconstrained. In general, each jet-parton assignment 
yields two different solutions for p u z and both are tried. This doubles the number of \ 2 minimizations performed for 
each Ijt event. For the dil channel the kinematics is actually under-constrained due to the presence of two neutrinos 
from each of the leptonic W-boson decays. In this case it is necessary to make an additional kinematic assumption 
in order to obtain a reconstructed invariant mass for each event. An integration over the undetermined additional 
kinematic variable (e.g. the azimuthal angles of the escaping neutrinos) is made and a reconstructed mass is calculated 
for each assumed value utilizing the kinematic fit. A weight is assigned to each of these possibilities by comparing 
the observed kinematics with those expected for the assumptions being made (e.g. the direction of the measured 
S}t)- The weighted distribution of reconstructed masses is then used to determined an m\ cco estimate for each event. 
Specific examples of this methodology are discussed in Sec. IIX.C.1I 

H. Event Selection Criteria 

After the full reconstruction software has been run, the following event selection criteria are employed to identify 
a sample of candidate it dil, Ijt, and had events. The three samples are statistically independent by construction. 
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3 In perturbative field theory the pole mass of a particle is defined as the pole of its renormalized propagator. 
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The selection of dil events begins by requiring a pair of high energy leptons, e + e~, , or e ± /j, T , each with Et 

exceeding a value around 20 GeV and |ry| < 2(1) for electrons (muons). The events are collected using the inclusive 
lepton triggers described in Sec. IV. Al and reconstructed using the algorithms described above. Several lepton selection 
criteria categories are usually employed ranging from "tight" , which yield a very pure sample of leptons, to "loose" , 
which can significantly increase the acceptance but with hadron mis-identification rates that are significantly higher. 
The analyses usually treat the various lepton-lepton categories (e.g. tight-tight, tight-loose, etc) separately since their 
signal-to-background ratio are significantly different. The leptons are required to be isolated. Backgrounds from Z 
events are removed by vetoing events for which either of the leptons, when combined with any other track in the 
event with Et > 10 GeV, yields an invariant mass consistent with the mass of the Z boson. The events are also 
required to have a large missing transverse energy, E^t > 20 — 25 GeV, and at least two jets with £j?* > 15 — 20 GeV 
(after correcting to the particle level) and \t)\ < 2.0 — 2.5. To reduce contributions from backgrounds for which the 
E^t arises from the mis-measurement of a jet or lepton energy, the angle between the E]t vector and any jet or lepton 
in the event is required to be larger than some minimum (e.g. A(j)(E^T — jet) > 20°). To remove backgrounds from 
b —> £cX — > ££qX decays the lepton-lepton invariant mass is required to be larger than 5 GeV/c 2 . If an e + e~ event is 
consistent with having originated from a photon conversion in material, or a fi + fj,~ event is consistent with a through 
going cosmic ray event, they are vetoed. Events with or 1 reconstructed jet are used as background dominated 
control regions to verify the accuracy of the background estimates. The total acceptance times branching fraction 
achieved is typically around 0.8% with a purity of around 70 — 75%. If a requirement that at least one of the jets is 
identified as a 6-jet the acceptance times efficiency falls to about 0.5% and the purity rises to about 90 — 95%. 

The selection of Ijt events begins by requiring a single high energy lepton, typically with Et > 20 GeV, collected 
using the inclusive electron and muon triggers described in Sec. IV. Al After employing the reconstruction algorithms 
described above, "tight" lepton selection criteria are used to identify the primary lepton in each event, which is 
required to be isolated. Events in which another isolated lepton is identified with Et > 10 GeV are vetoed in 
order to remove dilepton backgrounds. Dedicated algorithms are used to remove e events consistent with having 
originated from photon conversion in the material. The events are also required to have large missing transverse 
energy, E^t > 20 — 25 GeV, and at least three jets with E^ 1 > 15 — 20 GeV (after correcting to the particle level) and 
\rj\ < 2.0—2.5. Events with 1 or 2 reconstructed jets are used as background dominated control regions to verify the 
accuracy of background estimates. The total acceptance times branching fraction achieved is typically around 7 — 8% 
with a purity of around 15 — 20%. To improve the purity, most it analyses using the Ijt final state require that at 
least one of the jets in the event be identified as a 6-jet. Doing so yields an acceptance times branching fraction of 
around 4% with a purity of about 60%. 

The selection of had events begins with a data sample collected using the multi-jet trigger described in Sec. IV. Al 
After utilizing the reconstruction algorithms described above, events with an identified high energy lepton are removed. 
Since no missing transverse energy is expected in the all-hadronic final state, the observed E}t is required to be small 
relative to the uncertainty with which it is measured, E}t I Jv^t < 3. Finally, the events are required to have 6-8 
reconstructed jets, each with £j?* > 15 GeV (after correcting to the particle level) and \r\\ < 2. Events with 4-5 
jets are used as background dominated control regions in which to verify the accuracy of the background estimates. 
The total acceptance times branching fraction is about 13% with a purity of about 0.3%. To improve the purity, 
&-jet identification algorithms are employed, as well as multivariate discriminates that use kinematic and event shape 
information to suppress the background relative to it events. After these additional selection requirements are made, 
the total acceptance times branching fraction falls to approximately 3% and the purity rises to about 9%. 

All the analyses described below begin by using the samples obtained after employing these event selections. Often, 
depending on the analysis, additional criteria are used to further improve the tt purity of the sample. A variety of 
methodologies are employed to account for the remaining background contributions as described in the next section. 

VI. BACKGROUND PROCESSES 

After all selection criteria, background events still remain and must be accounted for in order to extract the it 
physics parameters of interest. A variety of background processes contribute to each it final state and fall into two 
basic categories: physics backgrounds and instrumental backgrounds. Physics backgrounds are those processes that 
share the same final state as the it signal sample we are aiming to isolate. For example, the pp — >• W + bbqq process, 
with the W decaying to ev or /iz/, is a physics background to the it — > Ivqq'bb sample. Instrumental backgrounds are 
those processes which mimic the it final state of interest due to an instrumental effect resulting in a mis-identification 
of some of the final state objects. For example, the pp — >■ W + bbqq process, with the W decaying to ev or \iv, is 
an instrumental background to the tt —¥ £^vbl~vb sample when one of the jets is mis-identified as a high energy 
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lepton. Although the mis-identification rates are typically very small, < 1%, the instrumental backgrounds can still 
significantly contribute to the final selected samples due to the very large production cross sections for the relevant 
QCD and W+jets processes. In general, the acceptances for physics backgrounds are estimated using Monte Carlo 
simulations, while the instrumental backgrounds are estimated using data control samples. For background processes 
making small contributions to the final sample, or with well understood theoretical cross sections, the normalization 
is taken from theory calculations. For all other cases the normalization is taken from the data using background 
dominated control samples and then extrapolated to the final signal sample. A brief review of the background 
processes that contribute to the various it final state event samples and how they are estimated is provided here. 

The generic QCD process pp — > jets has a production cross section that is about 9 orders of magnitude larger 
than the it production cross section. The jets produced pre-dominantly originate from wcis-quarks or gluons, although 
6-quark jets are produced in a few percent of these events. These generic QCD processes are the dominant background 
for the it — > qq'bqq'b final state with the instrumental background originating from the mis-identification of a udsg-jet 
as a 6-quark jet being the most important. For the it — > ivqq'bb final state, the QCD background is sometimes referred 
to as the "non-W" background and accounts for roughly 1/5 of the total background. It is an instrumental background 
and is a consequence of the mis-identification of a jet as an isolated high energy lepton and the mis-measurement of 
the S^t which makes it fall into the selected sample. If 6- jet identification is required, a further mis-identification of 
one of the udsg-jets is also necessary in order for the event to survive all selection criteria. All these mis-identification 
rates are taken from data control samples and applied successively in order to estimate the contribution of this 
background to the final selected sample. There are large uncertainties associated with this background which arise 
from the relatively low statistics of surviving events in the control samples once the identification criteria are applied 
(ie. because the mis-identification rates are quite small) and from systematic uncertainties accounting for kinematic 
differences between the control samples and the signal sample. For the it — > £ + vb£~vb final state this background is 
negligible since two jets would need to be mis-identified as an isolated high energy lepton and the S^t would have to 
be mis-measured. 

The VF+jets processes have production cross sections about three orders of magnitude greater than the it process 
and are an important background to the it — > Ivqq'bb and it — > t + vbi~vb final states. When the W decays leptonically, 
there is a genuine high energy lepton and genuine S^t in the event. For the Ijt final state, prior to making a 6- jet 
identification requirement, these processes appear as physics background. Accordingly their acceptance is estimated 
from Monte Carlo, however their normalization is taken from the data since there are large theoretical uncertainties 
associated with the predicted production cross-section. These uncertainties arise because complete calculations of the 
PT + Sjets and W+4jets cross sections, including heavy flavor contributions, are unavailable and current estimates rely 
on a mixture of partial calculations at lower orders and parton shower Monte Carlo models to extrapolate to larger 
jet multiplicities. When 6-jet identification is required, this background is separated into two pieces: the W+ light 
flavor (W + If) contributions, where the additional jets all originate from wds-quarks or gluons, and W+ heavy 
flavor (W + hf) contributions, where at least one of the additional jets originates from a c- or 6-quark. The W + If 
contributions is then an instrumental background since one of the udsg-jets would need to be mis-identified as a 
6-quark jet. As mentioned above the mis-identification rate of the 6-jet identification algorithm is taken from data 
control samples to estimate the contribution from this instrumental background, sometimes referred to as the "mis- 
tag" contribution. The W + hf contribution to the background uses Monte Carlo samples to estimate the fraction of 
surviving events which contain c- or 6-quark jets and the efficiency of the 6-jet identification algorithm to determined 
the number of events surviving the full selection requirements. There are large uncertainties associated with the 
heavy flavor fractions determined from the Monte Carlo, which are cross checked in data control samples. The 
VT+jets processes are the dominant backgrounds to the it — > ivqq'bb final states. For the it — > t + vb£~vb final states 
these processes are an instrumental background since one of the jets must be mis-identified as a high energy isolated 
lepton in order to survive the full selection criteria. This background is estimated from the data starting with the 
singlc-lcpton+i^T sample and then applying lepton mis-identification rates determined from data control samples to 
estimate the number of "lepton fakes" which populate the final dilcpton sample. 

A related background, Z/7+jets, has a production cross section roughly a factor of ten smaller than the W+jets 
background. It contributes to the it — > ivqq'bb sample only if the Z/7 decays to e + e~ or fi + p.~ and if one of the 
leptons escapes undetected to give rise to fake S^t- This latter effect is dominated by the limited geometric acceptance 
of the detector and is thus estimated using Monte Carlo simulation and found to be a negligible contribution to any of 
the samples. For the Ijt sample the only significant contributions arises from Z/7 decays to t + t~ with the subsequent 
decay of one of the r leptons to t ^ ev T v e or t — »■ pv T v^ in order to generate a lepton and E}t to survive the selection 
criteria. Once the relevant branching fractions are included, this turns out to be a small background and is estimated 
from the Monte Carlo. If 6-jet identification is required, the relevant efficiency and mis-identification rates are taken 
from data control samples and applied to the surviving Monte Carlo events. For the dil sample this is potentially an 
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important background since Z/'y decays to e + e~ or yield two high energy isolated leptons. However to survive 

the full selection criteria the B]t must be mis-measured. To further supress this background analyses raise the S}t 
requirement for events with a di-lepton mass consistent with mz ■ This is an instrumental background whose residual 
contribution to the final sample is estimated from the data using clean samples of identified Z — > e + e~ or to 
parameterize the fake fijx distribution and extrapolate to the ti signal region. There are relatively large uncertainties 
associated with this background arising associated with the extrapolation from the low jet multiplicity region which 
dominates the data control sample, to the >= 2 jets region of the signal region. 

The diboson processes pp — > VV (V — W, Z, or 7) have production cross sections within a factor of two of the it cross 
section. When the W and/or Z/'y bosons decay leptonically these processes make moderate (~ 10%) contributions 
to the Ijt samples and large contributions (up to ~ 30%) to the dil samples . These are physics backgrounds whose 
estimates are taken from Monte Carlo and normalized to theoretical predictions calculated to next-to-leading order. If 
&-jet identification is required, the relevant efficiency and mis-identification rates are taken from data control samples 
and applied to the Monte Carlo. 

The production of a single top quark via the electroweak interaction, pp — ¥ tq, has a production cross section about 
a factor of two smaller than the it cross section and makes a small contribution to the it — » ivqq'bb final state. It is 
a physics background and is estimated using Monte Carlo samples normalized to the theory predicted cross section. 



VII. SIMULATION SAMPLES 



A variety of Monte Carlo event generators are employed to estimate acceptances, efficiencies, and kinematics for 
it and background events. A summary of which generators are used to model each physics process is given here. In 
several cases a comparison of multiple generators is used to assess systematic uncertainties. For the analyses described 
in this review, the Monte Carlo generators are always combined with simulation software that models the detector 
response in detail. The resulting simulated samples are then treated just like the recorded pp collision data, using the 
same reconstruction software and particle identific ation algorithms descr ibed in Sec. El For all samples the decay of 
tau leptons is handled using the Ta uola package (Davidson et al , 2010() and the decay of c- and &-quark hadrons is 
handled using the Evtgen package ( Lange and Rvd . 19981 ). The parton shower and fragmentation processes required 
to simulate the evolut ion of a parent quark or gluon into a final state jet are modeled using dedicated sub-routines 



in either the Pythia (|Sj5strand et al ., 2001) or Herwig (Corcella et al 



2001) Monte Carlo software packages. The 



parton shower is modeled using DGLAP evolution (jGribov and Lipatovl . Il972l ) while the fragmentation process is 
simulated using a string model in Pythia and a cluster model in Herwig. Initial and final state QED radiation 
effects are included for all initial and final state charged particles. 

The CDF and D0 collaborations each employ several Monte Carlo generators to model it events. The Pythia and 
Herwig programs are general purpose event generators widely used throughout high energy physics. They model a 
wide variety of interactions and processes using leading order calculations (LO). Although the programs include only 
LO (e.g. pp — > it), the parton shower and QED radiation modeling can yield more complicated event topologies (e.g. 
pp — > tt+jets or ii + n f). However, the resulting cross sections and kinematics of these higher order final states usually 
carry large uncertainties since their modeling can be quite sensitive to arbitrary choices of control parameters offered 
by each generator. Reasonable choices of these parameters are made by making detailed comparisons between the data 
and Monte Carlo samples in dedicated studies using data samples of inclusive jet and minimum bias events. For most 
of all the it analyses described in this review, the higher order contributions affect the physics parameters of interest at 
a level that is negligible compared to the sensitivity of the present measurements. Consequently the Pythia program 
is often used to model the it events at the Tevatron. A more rigorous approach expected to y ield improve d cross 



section and ki nematic modeling of m ore complicated final states is offered by the Madevent ([Alwall et all 120071 ) 



and Alpgen (Man gano et al. 



2003) programs. These programs explicitly calculate the lowest order contributions 
pp — s- it + 1 parton) and then use Pythia or Herwig to model the parton shower and 
To model an inclusive sample, as observed in the data, requires the summing of several 

lj plus 



to a given process (e.g. 
fragmentation processes 

samples. For example, to model a it+ > jets sample requires an Alpgen simulation of it + Oj plus it 
ti+ > 2j samples. Since the parton shower can generate additional final state jets in each sample, care must be 
taken to avoid double-counting. Common choices to remove this overlap include those described in (jAlwall et all 
20081 ) . Several of the analyses described in this review employ a series of Alpgen+ Pythia or Ma devent+Pythia 
samples to simulate the inclusive it signal samples. The MC@NLO ([Frixione and Webberl . 120021 ) program is a full 
next-to-leading order generator used to simulate it samples employed for systematic studies. 

The W^+jets and Z/7+jets samples are generated using Alpgen with the parton shower and fragmentation modeling 
provided by Pythia. Sometimes the Pythia generator is used by itself to model the fully inclusive sample. In all 
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cases the normalization is taken from the data using selection criteria expected to yield a sample of events dominated 
by these processes. 

The diboson processes, VV (V = W, Z, or 7), are modeled either using ALPGEN+ Pythia or Pythia by itself. 
The samples are normalized to the theory predicted cross sections calculated to next-to-leading order. 

Electroweak pro cesses which result in the production of a single top quark are generated using either Singletop 
( Boos et all [2006) or Madevent with the parton shower and fragmentation modeling provided by Pythia. The 
samples are normalized using the next-to-leading order theory predicted cross section. The s-channel and t-channel 
processes are generated independently. 

Generic QCD processes are modeled using Pythia or Herwig or Alpgen+ Pythia. In general the predictions for 
the QCD background contributions to the final selected sample are taken from data control samples. However, several 
parameters important for tt analyses are taken from dedicated studies performed using multi-jet samples dominated 
by genuine QCD production processes. These studies often require some inputs from these QCD Monte Carlo samples. 
For example, the efficiency of the 6-jet identification algorithms is determined using a di-jet sample and the ratio of 
the data determined efficiency to the Monte Carlo predicted efficiency is used to correct the Monte Carlo predicted 
efficiency in tt samples. Another example, di-jet samples are used to compare the heavy-flavor contribution between 
data and Monte Carlo in order to derive a correction for the Monte Carlo predicted ratio W + hf/W + jets. The 
QCD samples used in these studies must be simulated in a manner consistent with that used to model the it and 
other background Monte Carlo samples so that the corrections derived can be used in the final selected samples with 
confidence. 



VIII. SOURCES OF SYSTEMATIC UNCERTAINTY 

A variety of systematic uncertainties are evaluated for each of the analyses described in this review. They fall into 
three general categories: uncertainties associated with the background modeling, uncertainties associated with the 
signal modeling, and uncertainties associated with the methodology employed to extract the physics parameter (s) 
of interest. The sources of uncertainty affecting the signal and background modeling and the methods employed in 
evaluating them, are common across all the it analyses and are summarized here. 

A. Systematic Uncertainties Associated with Background Processes 

The systematic uncertainties affecting the background modeling include instrumental effects related with the var- 
ious mis-identification rates required to promote a background physics process into the signal region, theory related 
uncertainties such as variations of the factorization and renormalization scales, k-factors related to NLO corrections, 
and differences in parton shower and fragmentation modeling, as well as uncertainties in normalization arising from 
uncertainties in lepton, 6-quark jet, etc. identification efficiencies. Several sources of uncertainty affect multiple 
background processes and their effects are accounted for in the analyses in a correlated manner. 

The mis-identification rate for hadrons faking an electron or muon signature and for udsg-jets faking a fc-jet are 
derived from data control samples. The mis-identification rates are parameterized by an empirically determined 
functional form, which includes the dependence of the rate on relevant kinematic variables (e.g. or \rj\). The 

validity of the parameterization is tested using statistically independent data control samples. Differences between 
the predicted rate, as determined using the parameterization, and the observed rate are used to assign systematic 
uncertainties on the mis-identification rates themselves. These uncertainties are then propagated to the number of 
predicted background events for those background processes for which the relevant mis-identification is required in 
order for that process to mimic a tt signal event. Since these mis-identification rates depend on the kinematics of 
the events, they also introduce a systematic uncertainty on the shape of the resulting kinematic distributions of 
background events surviving the full event selection criteria. These "shape" systematics are usually evaluated by 
re-weighting the relevant kinematic distributions using the mis-identification parameterizations varied within their 
uncertainties. 

The estimate for the Z/7+jets background in the dil channel is normalized using Z — >• tt events identified in 
the data. These events predominantly have small S^t and low jet multiplicity. The observed E}t distribution is 
parameterized and extrapolated into the tt signal region. The uncertainties on this extrapolation are dominated by 
the limited number of events with large S}t and the degree to which additional jets in the event affect this estimate. 
These uncertainties are added in quadrature and propagated to the relevant background estimate. Since the data 
events are selected requiring m« be consistent with mz, an additional correction is derived from MC samples to 
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estimate the number of Z/7 background events outside the mz region. Comparisons of different generators are used 
to assess a systematic uncertainty. 

The estimate for QCD background processes is data driven, taken from event samples collected using inclusive jet 
triggers. The QCD processes are a significant source of background in the Ijt and had channels. As discussed in Scc. lVIl 
they are instrumental backgrounds since they require the mis-identification of multiple objects in order to survive the 
relevant selection criteria. The rejection factors for each source of mis-identification are usually assessed individually 
and their product is taken to estimate the QCD background yield. This methodology assumes the rejection factors are 
uncorrelated with each other. Studies are performed to test the validity of these assumptions and assign systematics 
uncertainties. 

When generating the W+jets samples there is an ambiguity associated with the choice of factorization and renor- 
malization scales, represented in the literature as /xj?, /i^,, and/or Q 2 , at which to evaluate the relevant matrix element. 
Some typical choices for VF+jets processes are Q 2 = Pxwi ( m w +Pt w)i or ( m w + J2i Pt »)j where the sum runs over 
all the outgoing partons at a given vertex. The nominal background estimates are estimated using the last choice. 
A full set of additional Monte Carlo samples are generated using other alternative choices of Q 2 and the background 
predictions and kinematic distributions are re-evaluated. Differences are assigned as a systematic uncertainty affecting 
the W^+jets background processes. It should be noted that this ambiguity affects all the generated samples. However 
in the other samples it produces negligibly small systematic uncertainties. 

All the background processes are generated with at least two different Monte Carlo generators. For those Monte 
Carlo samples which reasonably model the observed data distributions, the background yields and kinematic distribu- 
tions are evaluated in full. Differences between the nominal Monte Carlo sample and any viable alternative sample are 
assigned as systematic uncertainties. This systematic uncertainty primarily accounts for differences in parton shower 
and fragmentation modeling. 

Most of the Monte Carlo samples employed in the analyses described here are based on leading order theory 
calculations. The inclusion of higher order corrections affect the predicted production cross section for the various 
background processes. These effects are accounted for in the analysis either by normalizing the Monte Carlo yields to 
next-to-leading order calculations, or by normalizing to the data itself. Some residual higher-order effects still remain. 
Most importantly is the W + hf fraction of the W-^+jets background. The fraction is taken from the Alpgen generator, 
but scaled to the heavy flavor fraction measured in the W + 1 jet sample. There is a large uncertainty associated 
with extrapolating that measured fraction to the W+ > 2 jets regions used in the tt analyses. The uncertainty is 
evaluated by comparing the measured heavy flavor fraction derived using different methodologies as well as comparing 
the result derived using the W + 1 jet sample to that derived using the W + 2 jet for the Ijt channel. The differences 
are assigned as systematic uncertainties and propagated to the predicted number of W + hf events contributing to 
the final background. 

The electron and muon identification efficiencies, the relevant trigger efficiencies, and the 6-quark jet identification 
efficiencies all have associated statistical and systematic uncertainties. These uncertainties are propagated through 
each of the affected background processes by varying each in turn within their uncertainties. These variations result 
in uncertainties on the number of contributing events for the affected background processes. For those efficiencies that 
have a kinematic dependence, their effect on resulting kinematic distributions is also evaluated by re- weighting the 
affected distributions using parameterizations that bracket the uncertainties of the relevant kinematic dependencies. 

The uncertainty on the jet energy scale is propagated by varying the parameterized jet energy corrections within its 
uncertainty and re-evaluating the predicted number of background events and the resulting kinematic distributions. 
This uncertainty affects all the background processes. 

The total uncertainty on the predicted background events is the quadrature sum of these various sources. When 
a particular source of uncertainty affects multiple background sources, the resulting correlation is accounted for by 
adding the resulting uncertainties linearly and then adding the resultant sub-total quadratically with the rest of the 
systematic uncertainties. 

B. Systematic Uncertainties Associated with tt Modeling 

There are systematic uncertainties affecting the modeling of it production and decay: variations of initial and final 
state radiation, variations in the parton distribution functions, variations in the parton shower and fragmentation 
modeling, variations in the modeling of color reconnection effects in final state interactions between the t and t decay 
products, and uncertainties associated with the trigger, lepton identification, and 6-jet identification efficiencies. Some 
of these sources of uncertainty also affect background processes and their effects are accounted for in the analyses in 
a correlated manner. 
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Uncertainties in the modeling of initial state radiation (ISR) and final state radiation (FSR) are evaluated by 
generating ti Monte Carlo samples with their ISR and/or FSR parameters varied within their uncertainties as de- 
termined using pp — > Z/j — » samples. Like it production at the Tevatron, this process is dominated by qq 

annihilation. The plj, M and the Nj et distributions are used to constrain the parameters in the Monte Carlo affecting 

the ISR modeling. The constraints are derived in bins of to?,, and extrapolated to to?,,, = 4m?. Since in the Monte 
<-> mm ± mm j 1 1 I — ■ — I 

Carlo generators the ISR and FSR processes are both modeled using DGLAP evolution (jGribov and Lipatovll 19720 . 
the same variations are used to assign systematics for FSR. Samples are produced varying ISR only, FSR only, and 
ISR/FSR simultaneously. The sample resulting in the largest difference relative to the nominal sample is used to 
assign the associated systematic uncertainty. 

The parton distribution functions r esult from multi-dimen sional fits to dozens of me a surem ents made by a variety 
of different experiments. The CTEQ (jPumplin et all 120021) and MRST ([Martin et all l2007t ) collaborations provide 
eigenvectors representing ±lcr variations of 20 uncorrelated parameters affecting the resulting PDFs. The difference 
associated with each variation is evaluated by reweighing the it events in the nominal Monte Carlo sample using 
each eigenvector in turn. The +ler and —la variations are quadratically summed separately. In addition, samples 
generated with different a s and different PDF fits are also generated and differences included in the associated PDF 
systematic uncertainty. 

The signal tt samples are generated with several different Monte Carlo generators. For those Monte Carlo samples 
which reasonably model the observed data distributions, the tt yields and kinematic distributions are evaluated in 
full. Differences between the Monte Carlo samples are assigned as a systematic uncertainty that primarily accounts 
for differences in parton shower and fragmentation modeling but also includes possible effects from differences in the 
modeling of the it pr spectrum, it spin correlations, and final state radiation. 

Since top quarks and W bosons decay quickly relative to the timescale associated with the parton shower and 
fragmentation processes (i.e. l/r t , 1/T\y << 1/Aqcd) it is possible that the products from the different top-quark 
decays could interact with each other via color reconnections. By default, all Monte Carlo generators ignore such 
reconnections since they are non-existent or irrelev ant for most processes. These color reconnection effects were first 
investigated at LEP2 (e.g. ( Abbiendi et all 2006f )). Since Tevatron's pp initial state carries the color charge, the 



situation and its modeling is significantly more complicated that at LEP, which used an e + e~ initial state. Recently, 
Monte Carlo models providing an adequate description of Tevatron data and including color reconnection effects have 
become available. The associated systematic uncertainty is evaluated by generating ti samples with these various 
color reconnection models enabled and comparing to the nominal Monte Carlo sample. Difference are assigned as 
systematic uncertainties. 

The uncertainties associated with the electron and muon identification efficiencies, the relevant trigger efficiencies, 
the 6-quark jet identification efficiencies, and the jet energy scale uncertainties are all evaluated in the manner described 
in the sub-section above. These uncertainties also affect some or all of the background processes. 

The total uncertainty on the predicted it yields and resulting kinematic shapes is the quadrature sum of these 
various sources. Correlations with the background processes are included. 



IX. MEASUREMENTS OF TOP-QUARK PROPERTIES 
A. Production Properties 

1. it Cross Section 

At hadron colliders top quarks are predominantly produced in ti pairs via the strong interaction. Measuring the 
ti production cross section is interesting for several reasons. First it allows a prec ise test of the predict ions from 
perturbative QCD. Computations of the pp — > it total rate are available at NLO ( Nadolskv et al\ . |2008). Several 
approximate NNLO computations, where the dominant (next to) next to leading logarithm terms are resummed to 
reduce sensitivity to the renormalization and factorizatio n scales, are also available ( Ahrens et al . 2010t Cacciari et al. 



2008 ; Kidonakis and Vogt . 20081 : lLangenfeld et al , 20091) . For a top-quark mass of 172 GeV/c 2 , the latest calculations 
predict a cross section of about <r{pp — > it) = 7.6 pb with a relative precision of about ±8 %. Comparing the 
measurements with these predictions is important since various BSM theories postulate new particles that couple 
preferentially to the top quark and thus predict a higher it productio n rate. A typical example is the production of 



new resonant particles that decay into it (jFrederix and Maltoni . 

It is also particularly important to measure the production cross section in different top-quark decay channels since 
new physics contributions can affect the various it final states, differently. Examples of such models are Two Higgs 
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Doublet Models among which are SUSY models ( Chung et al. , 2005 ; iGrossman , 1994 ) predicting the existence of 
charged Higgs bosons, H ± . If such a charged Higgs boson is light enough (mjj+ < m t — rrib), the decay t —> H + b can 
compete with the SM decay t — > Wb; the measured it production rate will then differ among the various final states 
due to contributions from charged Higgs decays. 

To make precise measurements of the it production cross section requires a good understanding of the reconstruction 
and identification efficiencies as well as a careful evaluation of the background processes that mimic the it signal. 
Because of this, the cross section analysis also then serves as a foundation on which most the rest of the it analyses 
are built. The it cross section measurement is performed by evaluating the following formula: 



N -N b 
e tI BjCdt 



(5) 



where N a is the observed number of events after selection, jVj is the estimated mean number of background events, e t j 
is the signal efficiency evaluated using it Monte Carlo, B is the relevant final-state-dependent it branching fraction, 
and J Cdt is the integrated luminosity for the particular set of triggers used for the measurement. It is worth 
noting that since e t j is increasing with the top-quark mass the measured cross section is quoted at a given m t . The 
numerator in Eq. [5] can be evaluated either using event counting after applying the final selection criteria or by fitting 
a discriminant variable that separates signal and background to estimate the relative contributions of each to the 
selected event sample. The use of 6-tagging is helpful to discriminate it signal from background. If 6-tagging is 
employed, N Q — Nb is usually evaluated by event counting since events with 6-jets have higher signal purity. 



tt Cross Section in the Lepton+jets Channel 

In the Ijt channel one of the W bosons decays to yield a high energy electron or muon plus a neutrino, while the 
other decays into two quarks. The final state thus contains at least four jets in the event, two of which are 6-jets, plus 
a high energy electron or muon and large S^p- This final state has the advantage of a reasonable branching fraction 
while achieving a good signal-to-background ratio. The event selection begins as described in Sec. IV. HI For the cross 
section analyses additional selection criteria are often applied. For example, the total transverse energy of the event, 
Ht — ^2 1 Ejp + X^jct -^t* + ^< T ' wnere the sums are over all high energy leptons and jets, respectively, is required 
to be above some minimum. At this point, further enrichment of the selected sample in it events is achieved by 
applying either a multivariate discriminant using topological information or using 6-jet identification. The topological 
information in the multivariate discriminate usually excludes any information explicitly related to 6-jet identification 
and thus yields a cross section measurement that makes no assumption about the flavor of the quark in the t — > Wq 
decay. Experimentally the final event samples are split according to lepton flavor (e or /i), jet multiplicity (3 or > 4 
jets) and (if used) number of identified 6-jets (0, 1 or > 2) and treated separately. Events with < 3 jets are used as 
background dominated control samples. 

The main background processes that can mimic a it event were discussed in Sec. IVI1 The most important contri- 
butions for the Ijt final state are, in decreasing order of importance, VU+jets production, QCD multi-jet processes, 
Z/j+jets diboson processes and electroweak single top-quark production. The Z/7+jets diboson, and single top 
backgrounds are evaluated using MC normalized to (N)NLO theoretical cross sections. As no reliable estimation of 
QCD multi-jet processes can be obtained by MC, this background is estimated using data. No reliable NLO cross 
section calculation exists to estimate the W+jets cross section at the large jet multiplicities relevant for it analyses. 
Therefore the normalization is derived from data and the MC used to estimate the VT+je ts kinematic d i stribut ions. 

The QCD background is estimated from data. For example, D0 uses a matrix method ( Abazov et al. , 2006bf) that 
consists of solving a 2D system of linear equations: 

Nloose — NqcD + N ther 

N tight = ZQCdNqcD + CotherNother (6) 

where Ni oose is the number of events in a data sample surviving loose lepton selection criteria, while N t i g ht is the 
number of events in a subsample surviving tight lepton isolation requirements. Nqcd is the number of QCD events in 
the loose sample and N other is the number of events that do not come from QCD events (dominated by W+jets and it 
events). The efficiency for the loose lepton from a lU-boson decay to survive the tight lepton selection criteria is e ther, 
which is estimated using W+jets and it MC corrected for any data-MC difference. Similarly, eqcd is the expected 
fraction of QCD events in the loose lepton sample that are expected to survive the tight lepton selection criteria. 
It is measured in a data sample with low S^t which is dominated by falsely isolated leptons after subtracting real 
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isolated lepton contributions. Solving this linear system yields an estimate for Nqcd- CDF uses a slightly different 
method to obtain an estimate of the QCD background contributions. A binned likelihood fit to the shape of the 
S^t distribution is used to determine the relative contributions of QCD and VF+jets processes. The S^t distribution 
for VF+jets processes is taken from the Alpgen MC, while the distribution for QCD processes is formed from data 
events that fail at least two lepton identification criteria. The fraction of QCD events in the high Jfyp region is then 
extrapolated from a fit to the low B}t region. 

The total amount of W+jets contribution is normalized to the number of events in data minus the estimated 
number of it, QCD, and electroweak events. Among the W+jets events the contribution from W + hf is found to 
be underestimated in the MC simulation. An additional scaling determined from data is applied to the W + hf MC 
samples. The W + hf content is measured in data events with < 2 jets using a neural network (NN) trained to 
separate W + If from W + hf. By counting the number of events in the part of the NN distribution dominated 
by W + hf or performing a fit to the full NN distribution using templates for udsg-jets, c-jets, and fe-jets, the ratio 
W + hf /W + If is determined after unfolding the relevant identification efficiencies. This can be compared with the 
ratio predicted in MC and a correction applied to the MC if necessary. In practice the CDF and D0 experiments find 
that the data-to-MC ratio for this heavy flavor fraction varies between 1.0 — 1.6 depending on the method employed 
and the event samples used. Their analyses typically use a value of around 1.3 with systematic uncertainties assigned 
that cover the full spread. This correction factor is then applied to all the jet multiplicity bins. D0 is using an 
additional scale factor for the W + c samples based on next leading order computations. 

The topological cross section measurements use the output of a multivariate discriminant, a neural network (NN) 
for CDF or a boosted decision tree (BDT) for D0, to further discriminate it signal from the background. These 
discriminants exploit differences in the kinematic properties of signal and background. Indeed because of the large 
top-quark mass, it events are on average more energetic, more central, and more isotropic than VF+jets and QCD 
events. Hence several kinematic distributions that are well modeled by MC simulation and that have a good separation 
power between signal and background are used as input to the multivariate discriminant such as Ht, the aplanarity, 
the sphericity (defined below), the lepton- ^j-jets invariant mass, the minimum angle between two jets, the minimum 
invariant mass of any two jets, etc. Multivariate output templates for the it signal and for the QCD, VF+jets and 
diboson backgrounds are then formed and used in a likelihood fit of the observed NN or BDT data distribution to 
determine the fraction of it events in the select ed samp l e. Th is fraction can then be used to determine the it cross 
section. D0 uses a nuisance parameter method (|Sinervol . 120031) to simultaneously extract the result and constrain the 
corresponding systematic uncertainties. 

Utilizing only topological information and a data set of 4.6 fb~ 
7.71 ± 0.37 (stat) ± 0.36 (syst) ± 0.45 (lumi) pb while D0 uses 5.3 fb" 1 to measure (lAbazov et al. 



CDF measures (lAaltonen et al 

-l 



2010c) a tJ = 



2011dJ) a fJ = 



7.68_ ; 64 (stat + syst) pb, both for a 172.5 GeV/c 2 top-quark mass. These measurements are systematics limited. 
Apart from the uncertainty on the luminosity, the dominant sources of systematic uncertainty are the uncertainties on 
the jet energy scale, jet energy resolution, and signal modeling. The precision of these measurements is around 9 %. 
Figure [5] shows the distributions for the CDF NN and D0 BDT discriminants used in these topological analyses. 

Aside from using pure topological information, the data sample can alternatively be enriched in it events by 
identifying 6-jets using one of the methods described in Sec. IV.FI In that case, apart from the W + hf background 
estimated as described above, the contribution from W + If events when a udsg-jet is mis-identified as a &-jet has to 
be determined. Since this mis-identification rate is difficult to model using MC, it is determined using a data sample 
dominated by QCD events and parametrized as a function of the relevant kinematic variables (e.g. the jet Et and 
if). This probability is then applied to MC W + If samples to estimate the number of mis-identified W + If events 
contributing to the final 6-tagged sample. A maximum likelihood fit is performed on the data, using the predicted 
background as an input, to extract the it production cross section. Often times the likelihood includes nuisance 
parameters that allow, for example, the background normalization and 6-jet identification efficiency, to vary within 
their uncertainties. The likelihood is constructed to account for the correlations between the background and signal 
acceptances for each of the nuisance parameters. 

Utilizin g &-jet identific a tion al gorithms that exploit the long lifetime of 6-hadrons CDF uses 4.3 fb -1 of data to 



measure (lAaltonen et all l2010d) a-, = 7.22 ± 0.35 (stat) ± 0.56 (syst) ± 0.44 (lumi) pb while D0 uses 5.3 fb" 



to measure (jAbazov et all . l2011dl ) a t j — 8.13_ ; 90 (stat + syst) pb, both for a 172.5 GeV/c 2 top-quark mass. The 
dominant sources of systematic uncertainty come from the luminosity determination, the estimation of the W + hf 
background and the uncertainty on the efficiency for identifying b-j ets. D0 comb i nes the topological and fe-jet iden- 
tification methodologies to obtain a t j — 7.78^o!67 (stat + syst) pb (jAbazov et all l2011dl ). Using 6- jet identification 



algorithms that identify leptons from semi-leptonic 6-hadron decays, CDF measures cross sections consistent with 
those a bove but with larger uncertainties due to th e lower efficiency and purity achieved by these soft-lepton algo- 



rithms (jAaltonen et all . 12010ft lAcosta et all . l2005al) . Figure [5] shows distributions from the CDF and D0 analyses 
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that use lifetime information to identify 6-quark jets. 

The largest systematic uncertainty for both the topological and b-tag based measurements is the uncertainty on 
the luminosity. It originates from the uncertainty on the luminosity counter acceptance and from the uncertainty on 
the knowledge of the pp inelastic cross section and totals about 6%. It can be effectively removed by measuring the 
it cross section relative to a cross section whose theoretical predictions are known to much better than 6%. A good 
choice at the Tevatron is to normalize to the inclusive Z/j — » 11 cross section. These events can be identified using 
the same trigger as the it sample, offer reasonably large statistics (1-2 orders of magnitude more than the it sample), 
and can be reconstructed with high purity (> 95%). The expression for determining the it production cross section 
is then 



-y—t££ /theoretical ■ 



(7) 



measured 



Hence the luminosity uncertainty is traded with the theoretical uncertainty on the Z/j cross section. The inclusive 
Z/j — > It cross section is measured using consistent trigger requirements and lepton identification w ith the it cross 



section. The Z/~f — > 11 cross section times branching ratio measured by CDF (jAaltonen et al 



invariant mass range of 66-116 GeV/c is uz/f- 
a t j from ([7]) a theoretical cross section of (crz/y->tt) 



2010d) in the dilepton 
247.8 ± 0.8 (stat) ± 4.4 (syst) ± 14.6 (lumi) pb To extract 



theoretical 



= 251.3 ± 5.0 pb (jAbulencia et all l2007al) is used and 



all rel evant systematic uncertainty correlations between the tt and Zj^j cross sections are included (jAaltonen et al 



2010c! ) to yield a tI = 7.82 ± 0.38 (stat) ± 0.37 (syst) ± 0.15 (theory) pb for the CDF topological selection and 



a tJ = 7.32 ± 0.36 (stat) ± 0.59 (syst) ± 0.14 (theory) pb for the CDF 6-tagging selection, both used a 4.6 fb data se t 
and assumed a 172.5 GeV/c 2 top-quark mass. These two measurements are combined to give (jAaltonen et aZj . l2010a) : 
a t j = 7.70 ± 0.52 pb for a 172.5 GeV/c 2 top-quar k mass. The precision of the combined result is now 6.8%, which 
is mo re precise than the latest theory predictions ( Ahrens et al\ . 2010t Kidonakis and Vogt . 2008 ; Langenfeld et al. 
2009h . 



tt Cross Section in the Dilepton Channel 

In the dil channel both W bosons decay to leptons. The event then contains two high Et isolated leptons, large 
EiT ■> and two 6-jets. This channel is often sub-divided further according to the explicit lepton identification criteria 
employed. When both leptons are reconstructed as an electron or muon the sample is dominated by t — ¥ Wb — ¥ evb and 
t —> Wb —> [ivb decays {it channel). The acceptance can be significantly increased, particularly for Wb — > tv — >• hvv 
decays, by applying very loose criteria for one of the lepton legs (^track). Measurements have also been performed 
in the final state with a r lepton decaying into hadrons explicitly identified from the decay of one W boson and 
an accompanying electron or muon from the other VF-boson decay {It channel). In general the dil channel has the 
advantage of a good signal to background ratio even without using 6-tagging but suffers from a smaller branching 
fraction than the Ijt or had channels. 

In the dil channel, the main source of background comes from the production of electroweak bosons that decay to 
charged leptons. For the il channel, it arises from Drell-Yan processes, Z/7 — > £ + £~ , and diboson processes when the 
bosons decays lead to at least two leptons in the final state. These backgrounds are reduced by requiring > 2 jets 
and large E^t- In the £t and ^track channels, the main background comes from processes with jets mis-identified as 
an electron, r, or isolated track, or from muons from semi-leptonic &-quark decay. This occurs mainly in VF+jets or 
QCD events. The diboson backgrounds are evaluated using MC normalized to (N)NLO theory cross sections. The 
contribution of Z/7+jets backgrounds to dil events in which the leptons are of different flavors (e.g. e/i events) is also 
evaluated using MC normalized to NLO theory cross sections. The contribution of Z/7+jets to dil events in which 
the two leptons are the same flavor (ee or /x/z), VF+jets normalization, and the QCD backgrounds are evaluated using 
data control samples. The jet-to-lepton fake rates are computed in a background dominated sample orthogonal to 
the signal sample (QCD di-jets, 7 + jets, same sign dilepton samples). By applying 6-jet identification a very pure it 
sample can be identified. For these analyses the background after 6-tagging is estimated in the same manner employed 
for the Ijt channel. 

As it benefits from a favorable signal to background ratio, the it dilepton selection, at least in the case of two 
well identified leptons, relies on a small number of simple requirements. Selecting it events decaying to dileptons 
requires first an inclusive lepton trigger as described in Sec. |Vj two high Et isolated electrons or muons (typically 
with Et > 20 GeV) or one isolated electron or muon and o ne isolated h igh pt track and generally at least two jets with 



Et > 20 GeV. In the tr channel, a set of neural networks (jGaleal . l2007j ) developed to separate hadronically decaying r 



leptons and jets is used to identify a r lepton in the event. Large S^t is further required (typically S^t > 25 to 35 GeV). 
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The S!t threshold is raised if the dilepton invariant mass is in the range of mz- Other topological cuts like a cut on 
Ht sometimes replace the S^t requirement particularly for the ey, channel since it suffers from the Drell-Yan process 
only through double leptonic r decays of Z/7 — > t+t~. In the It and £track channels, to reduce the contribution 
from events with fake SJt due to misconstruction, requirements are made on the lepton-i^ and jet-i^r opening 
angle. In the £t and ^track channels, additional requirements are made to further suppress the QCD background and 
&-jet identification is often employed. In the 11 channel, the signal-to-background ratio after all selection criteria is 
about 3:1 prior to and 15 : 1 after requiring at least one 6-jet in the event. The analyses typically use simple event 
counting and minimize a likelihood function that includes nuisance parameters constraining background contributions 
and other systematic uncertainties to determine o t -. 

In the ftrack channel, using 1.1 fb _1 of data and for m t — 175 GeV/c 2 CDF measures ( Aaltonen et al. , 2009a) 



a t - t = 8.3 ±1.3 (stat) ± 0.8 (syst) ± 0.5 (lumi) pb without 6-tagging and a tJ = 1Q.5_ 13 (stat)l^ (syst) ± 0.6 (lumi) pb 
using 6-tagging. Combining these two measurements gives a t j = 9.6 ± 1.2 (stat)^^ (syst) ± 0.6 (lumi) pb. In this 
channel the main systematic uncertainties come from the uncertainty in the background determination for the result 
without 6-tagging and from the uncertainty on the 6-jet identification uncertainty in the second case. 



Using lfb" 1 of data, D0 measures in the 11 channel without 6-tagging for a 170 GeV/c 2 top-quark mass (jAbazov et al 



2009J) tr t j = 7.5+J; 2 (stat)^^ (syst)!^ (lumi) pb and in the It channel using 6-tagging and treating tt events con- 



taining a r lepton that is mimicked by a jet as background a t j — 7. 6+43 (stat)^'^ (syst)^Q g (lumi) pb. In the ££ 
channel, the main sources of systematic uncertainties come from the signal modeling and the jet energy scale while 
in the It channel, these come from the signal modeling and the background estimation. Figure [7] shows the predicted 
and observed BJt distribution in both CDF and D0 from these tt dilepton cross section measurements. 

Both collaborations have provided preliminary updates of the tt cross section in the dil channel using more than 
5 fb~ . The latest measurement from CDF in the tl channel for a 172.5 GeV/c 2 top-quark mass using 5.1 fb _1 
of data is a tJ = 7 AO ± 0.58 (stat) ± 0.63 (syst) ± 0.45 (lumi) pb without 6-tagging and a tl = 7.25 ± 0.66 (stat) ± 
0.47 (syst) ± 0.44 (lumi) pb using 6-tagging. The latest measurement from D0 in the 11 channel without 6-tagging 
for a 172.5 GeV/c 2 top-quark mass using 5.3 fb -1 of data is a t j = 8.4 ± 0.5 (stat)!^ (syst) j^'g (lumi) pb. Both of 
these measurements are now limited by the systematic uncertainties. 

tt Cross Section in the All Hadronic Channel 

In the had channel, both W bosons decay hadronically into two quarks. The final state consists then of at least 
six jets among which two of them are &-quark jets. This channel has the advantage of a large branching fraction and 
no undetectable final state particles but is challenged by a QCD multi-jet background that dominates the signal by 
several orders of magnitude in cross section. The use of 6-tagging in this channel is therefore essential. The purity of 
the sample can be further increased by using kinematical and topological characteristics of the tt events. 

The trigger and event selection is described in Sec.|V] After additionally requiring that two jets are identified as 6-jets 
the signal-to-background ratio is around 1:2. To further separate the all hadronic tt signal from the QCD background 
multivariate approaches are used. CDF combines several kinematical jet or event variables into a neural network 
while D0 uses a relative likelihood method. The choice of the input variables to these multivariate discriminants 
is based on the signal-to-background separation power and agreement between the data and the background MC 
model as evaluated in a background dominated control region. The input variables mainly measure the amount of 
energy in the event, and the topology of the event. For example, the event centrality, aplanarity, and sphericity are 
used. The centrality is defined as the scalar sum of the jet t divided by the sum o f the jet energies. The even t 
aplanarity and sphericity are defined through the normalized momentum tensor Ai (jBrandt and Dahmen . Il979h : 



Mij = ^aPfPj 1 /Sa \fi*\ 2 where p** is the momentum vector of a reconstructed jet a and i, j = 1,2,3 denote the 
three spatial components of the jet momenta. The three eigenvalues of this tensor are noted in decreasing order 
Ai, A 2 and A3 with Ai + A 2 + A 3 = 1. The aplanarity is defined as A = |A 3 and reflects the isotropy of the event. 
The sphericity is defined as S — |(A 2 + A3) . Furthermore di-jet or tri-jet invariant masses are also used. Additional 
variables like the 77 and (j) moments of a jet ([Aaltonen et al\ . 2010c) can be added as they show a good discrimination 



between quark-initiated and gluon-initiated jets, the former being the dominant case for tt events while the latter is 
the dominant case for QCD processes. 

The QCD background estimate in the had channel is taken from the data using events from the lower jet multiplicity 
bins. To measure the amount of background remaining after 6-tagging, CDF first evaluates the probability of tagging 
a udsg-jet in a sample with exactly four jets, which is QCD dominated. This rate is parametrized in terms of jet 
Et, the number of tracks associated to the jet, and the number of primary vertices reconstructed in the event. This 
tag rate is then used to predict the number of tagged events in the signal region taking into account the possible 
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correlation among the jets. The accuracy of this background model is validated by comparing the data with the 
estimated background in an intermediate NN region (0.75 < NN < 0.85) correcting for the small ti contamination 
that lies in this region. D0 creates a background sample from data by attaching \ow-Et jets selected from events 
with six or more jets to events with four or five jets ensuring a compatible phase-space configuration between the two 
sources of jets. This is achieved by requiring some amount of S^t in the four or five jets events with a low I^t/Ht 
ratio to reduce effects from mis-reconstruction or from events containing hard neutrinos. This model of background 
is validated by comparing the five jets events in data with the constructed four+one jet background model. The 
background in the final sample used for the cross section measurement is constructed with four+two jet events mixed 
with five+one jet events. The di fference between these two constructions is taken as systematic uncertainty. 

2010eh cr tJ = 7.2±0.5 (stat)±1.0 (syst)±0.4 (lumi) pb while D0 uses 
6.9±1.3 (stat)±1.4 (syst)±0.4 (lumi) pb, both for m t = 175 GeV/c 2 . 



Using 2.9 fb" 



1 fb to measure (Abazov et 



CDF measures (lAaltonen et 
2010 ft a~ 



The uncertainty on the background modeling, on the jet energy scale, and on the 6-jct identification efficiency dominate 
the systematic uncertainty on these measurements. Figure [8] shows the distributions of the multivariate discriminants 
employed by CDF and D0 in their ti cross section measurements in the had channel. 

D0 also performed a ti cross section measurement in the r+jets final state with the r lepton decaying hadronically. 
The t lepton identification criteria are similar to these used in the It channel described above. A neural network 
is built to separate ti signal from the main background coming from QCD events. The cro ss section is extracted 
from a fit to the entire output discriminant distribution. Using 1 fb -1 of data, D0 measures ( Abazov et al 
a tJ = 6.9 ± 1.2 (stat)io;7 (syst) ± 0.4 (lumi) pb for a 170 GeV/c 2 top-quark mass. 



2010dD 



tt Cross Section Combination and Ratio of Cross Sections 



CDF and D0 measure the ti cross section in almost all possible decay channels assuming the SM branching fractions. 
A comparison of the results is sensitive to potential BSM contributions that affect the various decay channel in ratios 
differing from the SM, while a combination of the results yields an improved precision on the ti cross section. 

Within uncertainties all the measured ti cross sections in the different final states agree with each other and are 
thus combined. To simplify the combination all channels are constructed to be statistically independent. To compute 
the comb ined cross sect i on, D0 defines a joint likelihood function multiplying the Poisson probability for the different 
channels (| Abazov et aLl , l2009a) . Additional Poisson terms are added to constrain the QCD background in the relevant 
channels. Each systematic uncertainty is included in the likelihood function as a nuisance parameter. Correlations 
between channels are taken into account by u sing the same nuisanc e parameter. Combining the Ijt, the It and the 
£t channels measured using 1 fb" 1 , D0 finds (jAbazov et all l2009ch a a = 8.18±^;|f pb for m t = 170 GeV/c 2 . CDF 



performs a combination of its latest preliminary measurements, using up to 5.7 fb 1 of data in the dil, had , topological 
Ijt, and &-tag Ijt channels, by forming a best linear unbiased estimate (|Lvons et al. , 1988; V alassi . 2003 ) and taking 



into account the statistical and systematic correlations. For m t — 172.5 GeV/c 2 the prelim inary combination yield s 



7.46^Q gp pb. Both of these combinations are in agreement with the theoretical computation ( Langenfeld et al\ . 2009) 



By the end of the Tevatron run, a precision of about 6% could be achieved. 

Calculating the ratio of a t j measured in different final states or with different numbers of 6-tagged jets probes 
for the presence of non-SM decays of the top quark. Some BSM theories which might give rise to such effects 
were discussed in Sec. [Xj An example of such an approach is the measurement of the ratio of top-quark branching 



fractions, Rh 



B(t->Wb) 



, which can be expressed in terms of CKM matrix elements. D0 performed a simultaneous 



B(t^Wq) > 

determination of Rb and a t j (| Abazov et all l2008hl ). In addition to applying 6-tagging, a topological discriminant is 



used in the > 4 jets and 6-tag subsamples to further constrain the number of ti events. Using 0.9 fb" , the result 
is: Rb — 0.97lgQg (stat + syst) for a 175 GeV/c 2 top-quark mass in agreement with the SM expectation. The main 
uncertainty comes from the limited statistics while the systematic uncertainty is dominated by the uncertainty on the 

6-jet identification efficiency. 

Another example of BSM top-quark decays is t — >• H + b ( Chung et al. . 2005 : Guasch et al. . 1995f l. A large H + — > 



t + v branching fraction would result in a larger fraction of ti events showing up in the £t channel and fewer events 
in other channels than expected by the SM. On the other hand, a lepto-phobic charged Higgs would lead to fewer 
events in the 11 channel compared to Ijt. The ratios of the ti cross sections as measured in the various decay channels 
was determined by D0 after taking correlations into account. The ratio of the 11 to Ijt channel was measured 



as (|Abazov et all l2009d) R, 



jii/ljt 



0.86^o 1? while the ratio of the It determined cross section relative to all the 



others was measured as R a = 0.97Iq 29, both of which are compatible with the SM expectation. These 

were used to set limits on tau-onic and lepto-phobic charged Higgs decays as well as on the parameter space of the 
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MSSM (jAbazov et al. . 2009i ). CDF has also reinterpreted its cross section results to set limits on the charged Higgs 

2006bl) and assuming five possible Higgs decay modes: t — > W + b, t — > H + b — > rvb, 
-> i*66 and t — > H + b where H + — > W + h and h — > 66. Limits are set in the plane 
If no assumption is made on the charged Higgs decay, an upper limit on 



production (jAbulencia et al 
t -> H+b -> csb, i -> #+& 
(i? + , tan /?) for several MSSM scenarios 



B(t -> //"+&) of 0.91 is set at 95 % confidence level. 



tt Differential Cross Section 



The large statistics tt samples now available at the Tevatron enable a more thorough study of tt production. In 
addition to providing good statistics, the Ijt channel also offers good resolution on reconstructed quantities, and good 
sample purity. Thus, the most constraining measurements of top-quark properties and searches for BSM physics in 
the top-quark sample predominantly come from the Ijt channel. Both CDF and D0 have used fully reconstructed 
Ijt events to measure differential cross sections. Several BSM theories like technicolor, topcolor, or models with extra 
dimensions can distort the tt invar iant mass spectrum o r the transverse momentum distribution of the top quarks 
compared to the SM expectations ( Frederix and Maltoni . 20091 : Hill and Parkel 1994 ). Measurements of differential 
cross sections also test the perturbative QCD heavy-quark predictions. 



CDF measures the tt differential cross section with respect to the tt invariant mass, da/dm t j (jAaltonen et aLl . l2009d ) 



using the Ijt channel with 6-jet identification. The event selection and background estimate is described above. To 
reduce the systematic uncertainty introduced by the jet energy corrections, which affects the measured m t j spectrum, 
the jet energy scale is constrained in situ by comparing the measured invariant mass of the hadronically-decaying W 
boson to the world average VF-boson mass. The reconst ructed background subtrac t ed m f j distribution is corrected to 



the parton level using a regularized unfolding method ( Hocker and Kartvelishvili , 19961) . Apart from the statistical 
uncertainty, the largest systematic uncertainty arises from uncertainties in the parton distribution functions since the 
tail of the m t j distribution is very sensitive to the large- a; contributions to the PDF. The resulting differential cross 
section shows no evidence of BSM physics. 

D0 measures the inclusive differential cross section for 



da/dpip (jAbazov et al. 



PP 



tt + X as a function of the top-quark px, 
2010a) using the Ijt channel selected as described above. A kinematic fit, which con- 



strains the W-boson and top-quark masses to fixed values, is performed for all possible jet-parton assignments and 
the combination yielding the minimum fit chi-squared is used to associate the leptons and jets to the top quark. 
The reconstructed p T spectrum is corrected to the parton level. The uncertainties on the shape of the p T spectrum 
are dominated by statistics while the total systematic uncertainty integrated over pip is around 10%. The NLO or 
NNLO perturbative QCD calculations agree with the measured cross section in normalization and shape while the 
results from the Alpgen+Pythia and Pythia generators describe the shape of the data distribution but not its 
normalization. Figure [9] shows the measured differential it cross sections with respect to m t j or the top-quark pt 
measured by CDF and D0. 



2. Forward-Backward Charge Asymmetry 



At lowest order in QCD, the SM predicts that the top-quark pair production in pp interactions is charge sym- 
metric. However radiative corrections involving either virtual or real gluon emissions lead to a small charge asym- 
metry (jKuhn and Rodrieol Il998l Il999h . This asymmetry has its origin in radiative corrections to qq fusion from 
interference of the box diagram with the Born diagram while the contribution to tt production from gluon fusion 
remains charge symmetric. The real emission corrections tend to push the top quark backward, opposi te to the pro- 



ton b eam direction, while the virtual corrections push it forward, in the direct i on of the proton be am (Bowcn et 



20061) . This asymmetry is predicted to be around 5% ( Almeida et al. . 20081 Bowen et al. . 2006 : Dittmaier et al. 



20071 ) while the measured asymmetry depends strongly on the region of phase space being probed. Due to charge 
conjugation symmetry, the charge asymmetry can also be interpreted as a top- quark forward-backwa rd asymmetry. 



exchange (jCheung et al. 



extr a-dimension models (iDiouadi et al. 



2009) 



Some exte nsions of the SM (iFerrario and Rodrigol . 120091 ) . such as Z' — » tt decays (iMal 



cawi et all. 119961). t-channel W 



a new color octet resonance (jAntunano et all . 1200 

2009t ) predict different asymmetries. Both CDF ([Aaltonen et al. 



Ferrario and Rodrigo, 2008J) or 
2008 d|) and 



D0 (jAbazov et all . l2008d) measure the forward-backward charge asymmetry using Ijt events with 6-jet identification. 



The asymmetry Aft, is frame dependent and can be defined either in the pp or tt rest frame. The tt rest frame is more 
difficult to reconstruct but the asymmetry in the pp rest frame is predicted to be about 30 % smaller (jAntunano et all 
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20081 ). In the pp rest frame the asymmetry can be written as: 

N t (cos 6 > 0) - 7V t (cos6> < q) 
N t (cos9 > 0)+N t (cos9 < 0) 



APP 
A fb 



(8) 



where Nt(x) is the number of top quarks satisfying requirement x, cos 9 = —Qe ■ cos a p , a p is the polar angle between 
the top quark with the hadronic lU-boson decay and the proton beam, and Qi is the lepton charge from the leptonic 
decaying W boson. In the tt rest frame the asymmetry can be defined similarly except that 6* is the production angle 
of the top quark in the it rest frame. This angle is related to the rapidity, y, of the t and i in the pp frame by 



Ay = y t - yi = 2 tanh" 




(9) 



where s is the square of the center-of-mass energy. So the it rapidity difference in the pp rest frame can be used to 
measure the production angle 9* in the it rest frame. As both Ay and cos 9* have the same sign, the asymmetry can 
be written 



- t N t (Ay > 0) 
fb N t (Ay > 0) 



N t (Ay < 0) 



N t (Ay<0Y 



(10) 



CDF measures the asymmetry in both frames using fits to the cos 9 and to the Ay data distributions after background 
subtraction. The background fraction is constrained in these fits using the estimates from the cross section analysis 
described above. The tt events are reconstructed using a kinematic fit that constrains the VF-boson mass and the 
rrit — mj difference. The combination of jet-parton assignments which minimizes the fit x 2 is used in the analysis. 
The result s are corrected to the parton level using a matrix inversion technique. Using 1.9 fb -1 of data CDF 



measures ([Aaltonen et all l2008dl) Af b = 0.17 ± 0.07(stat) ± 0.04(syst) and A% = 0.24 ± 0.13(stat) ± 0.04(syst) 



fb 



These results agree within 2 standard deviations with the SM expectation. Using the same methodology CDF 
has a result using 5.3 fb -1 . The resulting forward backward asymmetry in the it rest frame is A l j h = 0.158 ± 
0.075 (stat + syst), consistent with the SM NLO prediction within two standard deviations. As shown in Figure ITUI 
they additionally observe that the asymmetry is rising as a function of the tt invariant mass, m t j, and determine 



Ay b (m t t > 450GeV/c 2 ) = 0.475 ±0.114 (stat + syst) (jAaltonen et al.l , l2011aj) . which is about 3.4 standard devia tions 



away from the NLO prediction of 0.088 ± 0.013 obtained using the MCFM program (jCampbell and Ellisl . 119991) . A 
preliminary CDF analysis using 5.1 fb -1 of data in th e dil channel measures a total asymmetry consistent with that 
observed in Ijt events and which also grows with m t j (jAaltonen et al , 2011(J) . 

D0 measures the asymmetry in the tt frame and performs a simultaneous fit for the asymmetry and the background 
fraction using the Ay distribution and a likelihood discriminant designed to separate tt from TU+jets processes. In an 
effort to remain independent of any potential BSM contributions, they do not correct for reconstruction effects but 
provide a parameterizat ion of the resolution. The D0 measured asymmetry corrected for background contributions 
from 0.9 fb" 1 of data is (|Abazov et all l2008d ) A l j h oha = 0.12 ± 0.08(stat) ± O.Ol(syst). D0 derives from this limits on 
the mass of a leptophobic Z' decaying to tt as shown Fig. [TO] The result has been rece ntly updated using th e same 
analysis techniques and 4.3fb~ 1 of data to obtain Af b obs = 0.08 ±0.04(stat)± 0.01 (syst) (|Abazov et aZ.I . l2011bh . Both 
these results are consistent within two standard deviations of the SM expectation, which is estimated using a NLO 
QCD MC generator, applying the event selection criteria, and determining the resulting asymmetry. 

These analyses will remain interesting as the data sets increase since the Tevatron's pp initial state is a CP eigenstate 
and offers increased sensitivity to tt A^ effects relative to the LHC, a proton-proton collider. 



3. Fraction of gluon-gluon Fusion Produced tt 



At the Tevatron, a variety of production mechanisms can yield a pair of heavy quarks in the final state. Distinguish- 
ing between the various initial state mechanisms is difficult since information is lost as the heavy quarks hadronizc 
and form bound states. However, as discussed in Sec. [TlJ top quarks decay before they hadronize and the initial 
state spin information survives to the final state. This can be used together with other discriminating variables to 
determine the fraction of tt pairs produced via gluon-gluon fusion. At the Tevatron, the SM predicts that gg — > tt 
accounts for about 15% of tt production, with the remainder coming from qq — > tt processes. This fraction varies by 
about ±5% (absolute) due to un certainties in the gluon PDFs. The measured qq — » tt contribution can be affected 



by BSM particles decaying to tt (jHill and Parkel Il994t lLane and Eichtenl . 119951 IZhang et all l2000h . 
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CDF uses the Ijt channel with &-jet identification to measure the relative fraction of gluon-gluon produced it pairs 
using two different methodologies. The first method exploits the fact that gluons are more likely to radiate additional 
gluons so that gluon initiated processes result in a higher track multiplicity on average than quark initiated processes. 
The observed low-pr track multiplicity distribution is fit to a contribution from background, from qq — » it, and 
from gg — > it processes. The distributions used in the fit are taken from MC samples calibrated using data samples 
known to be gluon- rich or gluon-depleted (VF+jets and QCD di-jet events). Corrections are made for the number of 
high energy jets in the event, for background contributions, and for acceptance effects. Using this first method CDF 
measures cr(gg — > ii)/ a(pp — > it) = . 07 ± .14 (stat) ± 0.07 (syst) which corresponds to an upper limit of 0.33 at 
95 % confidence level (|Aaltonen et aZ.I . l2008bh . 

The second method uses kinematic information to discriminate between the two production mechanisms. A kine- 
matic fit is used to determined the best jet-parton assignment to use. This allows a full reconstruction of the final 
state and of the top quarks and their kinematics. A variety of variables are combined using a NN, trained using 
it MC samples. The variab les include several top-quark production ang l es, an d several that are sensitive to it spin 
correlations as suggested in ( Mahlon and Parke , 1997 : Parke and Shadmlll996l ). A fit including background, gg — > it, 
and qq — > it contributions, is pe rformed to the observed NN distribution. Using this second method CDF employs 
a Feldman- Cousins presc ription (Feldman and Cousins! . 1998 ) to set an upper limit of a(gg — > ti)/a(pp —> it) < 0.61 



at 95 % confidence level (jAaltonen et aZ.I . l2009gl ). The dominant systematic uncertainties on this measurement come 



from the background shape and composition and from the difference between LO and NLO predic t ions. 

The two measurements are combined to yield <r(gg — > ti)/a(pp —> it) = 0.07+q'qy ( Aaltonen et aZ.I . l2009gl ) consistent 
with the SM expectation. 



B. Decay Properties 



1. W-Boson Polarization Fractions 



Since the top quark decays before forming a bound state, the t — > W + b process provides an opportunity to directly 
probe the tWb interaction vertex. In particular, measurements of the polarization fractions of the resulting W bosons 
can be used to test for BSM tWb interactions. In the SM the fraction of W + bosons from top-quark decay produced 
with longitudinal polarization is about /o = 0.70, with left-handed polarization is about /_ = 0.30, and with right- 
handed polarization is about / + = 0.0004. Contributions from BSM physics can give ri s e to anoma l ous co uplings 



at the tWb vertex that change these polarization fractions ( Aguilar-Saavcdr a et al. 
examples include the SM extensions discussed in ( Bernreuther et alL 2004b! :" Nam . 



2002) 



2007: iKane et aZ.I . I1992T ). Some 



Measurements of the W- 

boson polarization fractions can be used to set constrains on the vector and tensor form factors, ff 1 , , f^, and 
/2 . These terms appear in the general form of the Lagrangian describing the tWb vertex. In the SM only /-f has a 
non-zero value. 

The VF-boson polarization affects several kinematic variables, which can be used to measure (/o,/+). The most 
commonly used is 9* , the angle between the VF-boson momentum direction in the rest frame of the top quark and 
the direction of the down-type fermion decay product of the W boson (charged lepton or d— , or s-quark) in the rest 
frame of the W boson. The VF-boson differential decay rate can be written in terms of cos 9* as: 



1 dr 



Yd cos 9* 



*) 2 + -/ (1 - cos 2 9*) + -/+(! + cos0*) 5 



(11) 



Estimating cos 9* requires each it event be fully reconstructed using, for example, a kinematic fit. The charged lepton 
Et is also sometimes used since it is correlated to the W-boson polarization. It produces a less precise measurement, 
but is insensitive to uncertainties in the jet energy scale. Both CDF and D0 have recently used cos 9* to measure the 
PF-boson polarization fractions with 2.7 fb _1 and 1 fb _1 of data, respectively. 

CDF uses the matrix element method described in Sec. IIX.C.1I to determine the polarization fractions in the Ijt 
channel with &-jet identification. In this analysis the LO it differential cross section is expressed in terms of cos 9* and 
the fractions (fo, /+) while the top mass is fixed at 175 GeV/c 2 . The fit to the data is performed under three different 
scenarios, a) a model independent fit to both f and /+ simultaneously, b) a fit to determine / + with fo fixed to its SM 
value, and c) a fit to determine fo with / + fixed to zero (its SM value given the current level of precision). Corrections 
for background contributions and corrections for the v ariation of the event s election acceptance as a function of 
(fo, /+) are included. The model independent fit yields i) Aaltonen et al . 2010gl ) f = 0.88 ± 0.11 (stat) ± 0.06 (syst) 
and /+ = —0.15 ± 0.07 (stat) ± 0.06 (syst) with a correlation coefficient of —0.59. These measurements are statistics 
limited. Fixing fo to its SM value yields /+ = —0.01 ± 0.02 (stat) ± 0.05 (syst) while fixing / + to zero yields 
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fa = 0.70 ±0.07 (stat)±0.04(syst). The dominant systematic uncertainties arise from uncertainties in the background 
modeling and from uncertainties in ISR/FSR and parton shower modeling in it events. 

D0 uses a kinematic fit to fully reconstruct the event and calculate cosO* in the Ijt and dil channels. In the Ijt 
channel events are reconstructed using a kinematic fit with fixed lU-boson and top-quark masses. In the dil channel, 
given a top-quark mass, a four- fold ambiguity exists for each set of jet-parton assignments considered. In the analysis 
all real solutions are considered. The resulting reconstructed cos 9* distributions are compared with templates derived 
from it MC with different lU-boson polarization fractions and taking into account background and reconstruction 
effects. For a top-q uark mass of 172.5 GeV /c 2 the measured lU-boson polarization fractions determined from a model 
independent fit are (jAbazov et adl2008dl) f = 0.43±0.17 (stat)±0.10 (syst) and /+ = 0.12±0.09 (stat)±0.05 (syst). 
These measurements are statistics limited. The main systematic uncertainties come from the background or tt 
modeling and from the MC template statistics. 

D0 combines the measure d VF-boson polarization frac tions with the measurement of the electroweak single 



top production cross sect ion (jAbazov et all l2007al l2008bl ) to set limits on the fi, and f% BSM form fac- 



tors ( Abazov et aa l2009h). 



2. Decays via Flavor Changing Neutral Currents 



In the SM quarks can change flavor only via lU-boson interactions, so that at tree level decays of the type, t — > qZ 
or t — > <77, where q — u, c, are forbidden. Thus the rates for these decays are suppressed because they can only 
occur through quantum loops involving the W boson. Because the top quark is so massive, the quantum loop 
corrections arising from new particle loops are much larger than is the case for the light quarks. This makes the top 
quark a potentially great place to look for evidence of BSM Flavor Changin g Neutral Currents (FCNC) decays. For 
example, the SM predicts B(t -> qZ) « 0(1O~ 14 ) (jAguilar-Saavedral [20M) . while new physics models like 2HDM , 
supersymmetry, and topcolor-assisted technicolor predict branching fractions as large as 0(1O -4 ) ( Larios et a/.l . l2006h . 
An observation of a branching fraction significantly larger than predicted by the SM would be unambiguous evidence 
of BSM physics. 

CDF has searched for the FCNC decay t — s- qZ using tt events when the top quarks decay as tt — > qZqZ or 



tt 



qZbW and one of the Z bosons decays into a e + e 



or 



pair ( Aaltonen et al. . 2008i| ). while the other 



boson in the event decays hadronically. The experimental signature is thus very similar to the Ijt channel except 
that the final state neutrino is replaced by a charged lepton. The dominant background process is the production 
of Z/7+jets. Smaller background contributions come from the SM tt, WZ, and ZZ processes. The event sample is 
divided into sub-samples depending on whether or not any of the jets have been identified as 6-quark jets. A kinematic 
fit is employed similar to the one described in Sec. IV.GI except that one or both of mw constraints are replaced by 
mz constraints depending on which final state hypothesis is being tested. To separate the FCNC signal from the 
background, requirements are made on the transverse mass of the Z and the four leading jets, the of the four 
leading jets, and on the kinematic fit % 2 . A binned likelihood fit to the \ 2 distribution is simultaneously performed 
for the sub-samples with and without an identified 6-jet and the fraction of FCNC decays in the data is extracted. 
In the fit the normalization for the dominant Z/7+jets background is allowed to vary while the tt cross section is 
fixed to the observed it — > Ivqq'bb cross section. The systematic uncertai nties a f fectin g the signal and background 
shapes and normalizations are taken into account using the technique of (jReadl . 1999). The dominant systematic 
uncertainty arises from the jet energy corrections an d in uncertainties related to t he fraction of Z/7+jets events that 
have identified &-jets. A Feldman-Cousins procedure (jFeldman and Cousind . ll998h is employed to set the limit. Using 



1.9 fb , CDF data is consistent with the SM expectation and an upper limit of B(t —> qZ°) < 3.7% is established 



at 95 % confidence level. A similar analysis has been performed by D0 searching for it —> qZqZ or it 
with a leptonic decay of the Z or W bosons. Using 4.1 fb -1 in the trilcpton final state, D0 does not ob 
of anomalous coupling and set a limit of B(t — > qZ°) < 3.2% at 95% confidence level (jAbazov et al 



qZbW but 



searches can also be performed using the electroweak production of single top quarks (jAbazov et al. 





2011eh 




2010g 



Similar 



3. Decays to Charged Higgs Bosons 



Many theories beyond the SM propose an extension of the Higgs sector to multiple Higgs doublets (Grossman! . 



1994 ). The simplest extension of the SM Higgs sector is the two Higgs-doublet model (2HDM) (jGunion et aZ.1 . 11989) . 



In the 2HDMs, the Higgs sector consists of three neutral Higgs bosons (h, H, A) and two charged ones (H ± ). Several 
types of 2HDMs exist and differ in their strategies to avoid large flavor changing neutral current effects. In type I 
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2HDM, only one of the Higgs doublets couple to fermions. In type II, one doublet couples to the up-type quarks 
and neutrinos and the other to down-type quarks and charged leptons. The minimal supersymmetric standard model 
(MSSM) is an example of type II 2HDM. In models of type III, both doublets couple to fermions and FCNCs are 
suppressed through other mechanisms. An observation of a charged Higgs boson would be unambiguous proof of a 
BSM Higgs sector. 

At the Tevatron direct production of a single charged Higgs boson can occur by qq' annihilation. Searches for 
qq' — > H + — > tb have be en performed by D0 assuming a charged Higgs boson mass larger than the top-quark 



mass (jAbazov et ali . l2009ij ). This channel has the same signature as the electroweak production of single top quarks. 
Direct production of H + H~ thro ugh the weak interaction has a rather small expected cross section, of the order of 
0.1 pb ([Carena and Haberl . 120031 ) and so no such search was performed so far. For nijj+ < m t in 2HDM, the top 
quark can decay to H + b, hence competing with the SM decay to W + b. As seen in section HX.Al the ratio of it cross 
sections as measured in different decay channels is thus sensitive to BSM top-quark decays. The measured it cross 
sections have been used to set limits in the plan e {m H + , tan /?) where tan /3 is the ratio of vacuum expectation values 
of the two Higgs doublets ( Abazov et ali [20091 ). At small tan/3 (< 1), the top branching fraction into charged Higgs, 
B(t — > H + b), can be large. In that regime, the charged Higgs boson predominately decays to cs for low charged Higgs 
mass (< 130 GeV ) or to t*b — > Wbb for large mass. At larger tan/3 (> 15), B(t — > H + b) can also be large, but in 
that case the charged Higgs boson decays almost 100 % of the time to t + v. 

Using 2.2 fb _1 , CDF has performed a direct search for decays of top quarks into charged Higgs bosons assuming 
H + — > cs using the Ijt final state including 6-jet identification. The di-jet invariant mass distri bution, which is 



otherwise assumed to come from the M^-boson decay, is used to search for evidence of such decays (jAaltonen et 



2009kJ). The tt events are reconstructed using a kinematic fitter similar to the one described in Sec. IV.Gl except that 



the mw constraint for the hadronically decaying VF-boson candidate is dropped and templates are built for signal 
and background processes. The signal is modeled using H + MC where the H + decay is forced to cs and T H + = 
is assumed for 60 < m H + < 150 GeV/c 2 . A binned likelihood fit, including systematic effects, is used to exclude 
B(t — » H + b) larger than 0.1 to 0.3 at the 95% confidence level assuming B(H + — > cs) = 1. 



C. Intrinsic Properties 

1. Mass 



The top-quark mass, m t , is a free parameter in the SM and must be experimentally determined. A precision 
determination of m t is important since quantum loops including top quarks contribute large corrections to theory 
predictions for many precision electroweak observables. For example, the SM predicts a precise relationship between 
the W- and Z-boson masses, 



mz 



J = (1 -sin 2 IV) (1 + Ap) 



(12) 



where sm9w is the weak mixing angle and Ap is at tree level and non-zero once quantum loop corrections are 
included. The dominant quantum corrections are quadratically dependent on the top-quark mass and logarithmically 
dependent on the Higgs boson mass, Ap = /(m 2 , lnm#). The experimental program then consis ts of measuring 



mz, and m t as precisel y as possible in order to constrain mu - detailed discussions can be found in (jFlacher et al. 



, 2009: 



Nakamura et a/.l .l2010). The experiments at LEP and SLAC precisely determined mz, the Tevatron and LEP2 experi 
ments precisely determine mw , and the Tevatron experiments alone precisely determine m t . It is i mportant to measure 
mt in all the different top decay channels since BSM contributions can affect them differently (jKane and Mrennal . 
19961) . 

For a free particle, the physical mass is usually defined as the pole of its renormalized propagator. However this 
definition is ambiguous when dealing with colored particles like quarks. Thus the definit ion of the top-quark pole 
mass is intrinsically ambiguous on the order of hoc o due to non-perturbative QCD effects (jBeneke and Braunl . 1 1994 ; 



Bigi et all 11994 ISmith and Willenbrockl . 119971 ). For an unstable particle (like the top quark), the pole mass enters 



the description of the resonance throug h a Breit-Wigner fu nction. Other mass defin itions exist such as the MS -mass 



using the MS renormalization scheme ( Bardeen et al. , Il978l) or those discussed in (jHoang and Stewart! [2008) . The 



MS-mass is only sensitive to short distance QCD effects and is often used to describe the mass of the light quarks 
when the typical energy of the process is much larger than the quark mass in question. At the Tevatron and LHC, the 
top-quark mass will be directly measured through the reconstruction of its decay products and the measured mass is 
taken to be the pole mass. However, aside from the theoretical ambiguity mentioned above, additional ambiguities 
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are introduced since the measurements are calibrated using MC generators that include model dependent descriptions 
of the parton shower and hadronization processes. These affect the experimentally reconstructed top-quark mass in 
a model depen dent way and introdu ce a set of modeling related systematic uncertainties that are of the order Aqq d 
or even larger ( Buckley et al. . lioilh . 

Aside from having to isolate a relatively pure sample of it events, there are a few experimental challenges to making 
a precise determination of m t . The most sensitive analyses fully reconstruct the it kinematics using a kinematic fit 
or a matrix element technique. In order to do so it is necessary to assign each final state jet as having originated 
from a particular final state quark. As discussed in Sec. IV. Gl this gives rise to a combinatoric ambiguity that dilutes 
the resolution of the reconstructed it kinematics. The resolution with which the it kinematics can be reconstructed 
is further diluted by the addition of ISR/FSR jets, the loss of quark jets due to limited geometric acceptance, and 
jet algorithm merging/splitting effects. Moreover, large systematic uncertainties can arise from uncertainties in the 
jet energy corrections. Roughly, each 1% of uncertainty in the jet energy corrections gives a 1 GeV/c 2 uncertainty 
in the top-quark mass. This was a limiting systematic uncertainty for early m t measurements. With larger statistics 
samples the jet energy corrections can be calibrated in situ by constraining the corresponding jet-jet invariant mass to 
the world average mw in events containing &W — > qq' decay. Doing so significantly reduces the associated systematic 
uncertainty, which then scales with the statistics of the it sample itself. Most measurements in the Ijt and had 
channels now take advantage of this in situ calibration by performing simultaneous fits to mt and the jet energy scale. 
Other measurements seek to identify experimental observables that are only weakly dependent on the jet energy 
corrections. 

For all analyses described below it is necessary to calibrate the measured m t using MC pseudo-experiments (ie. 
ensemble testing) in order to quantify the influence of the simplifications made in the employed techniques. Each 
pseudo-experiment is built from a mix of it signal and background events in the proportions estimated from the 
relevant cross section analysis. The number of events contributed by each physics process is allowed to vary according 
to a Poisson distribution. Once a set of pseudo-events have been assembled, they are then treated just like the data 
and a determination of m t is made. Then, from an ensemble of pseudo-experiments, the relationship between the 
true top-quark mass (ie. the m t value used when generating the it MC sample) and the mean measured mt can be 
established and, if necessary, corrections can be derived - usually as a linear function of the measured value. This is 
done separately for each analysis technique. The ensemble of pseudo-experiments is also used to study the expected 
statistical uncertainty and the statistical behavior of each analysis technique. Once a particular analysis has been 
calibrated, its performance is checked using MC samples. 

Four methods are used at the Tevatron to measure the top - quark mass. They are described below. For recent 
detailed reviews on top-quark mass measurements, see (Fiedler, 201fi Wicke . 2010h . 



Template Method 

This traditional method begins by choosing an observable correlated with the top-quark mass. Distributions of 
this observable are then constructed using MC samples generated with varying m t as input. The data distribution is 
then compared with these MC templates using a maximum likelihood fit. The observable most correlated with m t is 
the reconstructed invariant mass of the it decay products, m l t cco , which is the most often used and is estimated using 
the kinematic fits discussed in Sec. IV.GI If a given analysis uses the hadronic 14^-boson decay to calibrate the jet 
energy corrections, an observable correlated with mw is used - usually the invariant mass of the associated jets, rrijj 
- and the maximum likelihood fit is extended to two dimensions. The template method is relatively simple and can 
be easily extended to several decay channels or sub-samples. The observable(s) used as estimators can be chosen to 
minimize sensitivity to specific systematic uncertainties. The statistical sensitivity of this method is sometimes worse 
than other methods discussed because it does not use the full event information, nor take advantage of event-by-event 
differences to weight more heavily those events with kinematics resulting in improved m\ cc ° resolution 



A d etailed description of using the template method to measure mt in the Ijt channel can be found in (jAbulencia et al 



2006c ). The most recent published result uses 1.9 fb 1 of CDF data and simultaneously fits the Ijt and dil channels 



using th e in situ W — > qq' decay s to constrain the jet energy corrections and reduce the associated systematic uncer- 
tainties ( Aaltonen et all l2009eT ). In the Ijt channel a kinematic fit is used to determine (m r t cco ,mjj) for each event. 



Since the it kinematics are under constrained in the dil channel, a neutrino weighting algorithm described below is 
used to determine m\ cc °. It is used together with Ht in this particular analysis. Two dimensional templates are 
constructed for each channel and a joint likelihood fit is performed to measure mt ■ Templates are built from simulated 
samples of it and background events using MC generated at disc rete values of ro«. The shape of the templates at 
these discrete points is described using kernel density estimators (Scott, 1992 ; Silverman . 19981) and then smoothed 
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and interpolated using the method of ( Loader , 19991) to enable an estimate of the template shapes for any arbitrary 
value of rnt- The data in each channel is further divided into sub-samples of varying signal purity and m\ cco resolution 
using the number of reconstructed jets and the number of identified &-jets in the events. A full set of templates are 
created for each sub-sample. The two dimensional data distributions are compared to the resulting templates using 
an unbinned maximum likelihood fit to determine nit, the best-fit jet energy scale correction, and their associated 
uncertainties. The resulting top-quark mass is m t — 171. 9± 1.7 (stat + JES) ± 1.1 (syst) GeV/c 2 for the simultaneous 
fit of the Ijt and dil samples. The dominant systematic uncertainties arise from residual uncertainties in the jet 
energy correctiontQ and from variations of the MC generator used to model it events. An example of one of the m\ ec ° 
distributions is shown in Fig.[TTJ CDF has also used the template method to determine mt from a sa mple of Ijt events 



in which the 6-jets arc identified using one of the soft lepton tagging algorithms described in Sec. IV.Fl ( Aaltonen et al. 



2009i). 



CDF also employs the template method to determine m t using observables with minimal sensitivity to the jet energy 
corrections ( Aaltonen et al. . I2009il . l2010hh . Since the transverse momenta of the top-quark decay products depend 
approximately linearly on m t , the decay length of the 6-hadron measured in the transverse plane, Lt, or the Et of the 
lepton from the TU-boson decays can be used to build templates. These observables are almost completely independent 
of the jet energy corrections. Small dependencies remain due to uncertainties in it and background acceptances that 
arise from changes in the event selection efficiencies when the jet energies are varied. The shape of the Lt distribution 
is calibrated using di-jet data samples enriched in bb. The lepton Et scales are calibrated as discussed in Sec. IV.CI 
The two-dimensional Lt and lepton Et distribution is used to determine mt and the resulting likelihood curve is 
shown in Fig. HU Using 1.9 ftT 1 CDF measures with these observables mt = 170.7 ± 6.3 (stat) ± 2.6 (syst) GeV/c 2 . 
The measurement will clearly benefit from more statistics and has a very low systematic uncertainty coming from 
jet energy calibration. However it suffers from relatively large systematic uncertainties arising from the Lt and 
Et calibration. The precision of these calibrations should improve as the statistics of the relevant control samples 
increase. 



A de tailed description of using the template method to measure m t in the dil channel can be found in (jAbulencia et 



l2006al ). Recent measurements use essentially the same methodologies on significantly larger data sets. As described 
in Sec. IV.Gl since the kinematics of the dil events is under constrained it is necessary to make an additional kinematic 
assumption in order to obtain a reconstructed invariant mass for each event. The CDF and D0 template analyses in 

the dil channel primarily differ in the additional kinematic assumption they choose to make. 

[2009ll 



CDF 's most recent template analyses in the dil channel use the neutrino (/) weighting method (jAaltonen et al. 



2010ml) . which integrates over the azimuthal angles for the two neutrinos in each event. A x 2 function is minimized to 
estimate mJ eco for each event. This \ 2 function includes a term that constrains the measured quantities within their 
uncertainties to the assumed tt kinematics. Templates of the m r t cco distributions are obtained from simulated samples 
of tt signal and background events generated at discrete values of m t . The shape of the templates is parameterized 
and interpolated to enable an estimate of the template shape at any arbitrary mt value. The top-quark mass is 
determined using an unbinned likelihood fit to the data. In one of these analyses ( Aaltonen et aL 2010ml) the fit 



to the mj™ distribution i s extended to include information from the tot2 variable suggested in (jBarr et all . 12003 ; 
Lester and Summers! . 1999). CDF also develop ed a template based me asurement where the m t dependent SM tt cross 



section is included as an additional constraint ( Aaltonen et al .. 2008a). 



The most recent D0 template analyses in the dil channel use a neutrino weighting method and the so called matrix 



weighting method ( Abazov et al. . 2009el ). The neutrino weighting algorithm employed integrates over the neutrino 
rapidities as the assumed input to the kinematic fit. Weights are assigned by comparing the resulting neutrino 
momenta solutions to the measured S^t- The first two moments of the weighted distributions are used from each event 
to build m\ eco templates and extract from the data distribution the most probable top-quark mass. The templates 
are constructed in two different manners, using binned probability density histograms and using probability density 
fit functions, which are later combined. The matrix weighting technique integrates over assumed top-quark masses 
and solves for the t and t momenta. Each solution is then weighted by the probability to measure the observed lepton 
en ergy in the top-quark rest f rame given the assumed top-quark mass using the matrix element based expressions 
in (jDalitz and Goldsteinl . ll992l ). For each event m\ cco is taken from the maximum of the resulting weighted distribution 
and used to build templates. 

Usin g 3.4 fb -1 , CDF uses a combination of the mJ cco and mT2 distributions in the dil channel to measure (jAaltonen et 
12010ml) : m t = 169.3 ± 2.7 (stat) ± 3.2 (syst) GeV/c 2 . Using 1 fb~\ D0 combines the results from the neu- 



trino weighting and the matrix weighting techniques in the dil channel to measure (jAbazov et all l2009et ) m* 



4 The measurement treats the jet energy scale as a constant. The residual jet energy scale systematic accounts for the known variations 
of these corrections as a function of jet Et and 77. 
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174. 7± 4.4 (stat) ± 2.0 (syst) GeV/c 2 . These measurements are still limited by statistics. The systematic uncertainty 
is completely dominated by the uncertainty in the jet energy corrections as no in situ jet energy scale calibration is 
possible in the dil channel. 



A d etailed description of using the template method to measure m t in the had channel can be found in ([Aaltonen et al. 



2010eh in which CDF performs a simultaneous measurement of the tt cross section and the top-quark mass using a 



NN based event selection and requiring at least one identified &-jet. A kinematic fit is used and all possible jet-parton 
assignments consistent with the 6-jet information are considered. The combination yielding the smallest fit \ 2 is used 
to determine m\ cco , + , and m^~ for each event. Templates are built from the m\ cco distribution and from the 
nijj distribution for a set of MC samples generated at discrete mt and jet energy correction values. These templates 
are fit to a Gamma plus Gaussian function, and the resulting parameters are smoothed and interpolated to enable 
an estimate of the template shape at any arbitrary mt value. The top-quark mass is determined from an unbinned 
likelihood fit to the data that includes an in situ constraint on the jet energy corrections using the mjj distributions. 
Using 2.9 fb _1 , CDF measures m t = 174.8 ± 2.4 (stat + JES)]^ (syst) GeV/c 2 . The measurement is limited by 
statistics and the dominant source of systematic uncertainty arises from residual jet energy correction uncertainties 
and in method related systematics. 



Matrix Element Method 

The matrix element method is a more sophisticated method for measuring a particle property and constructs a 
likelihood curve for each event by comparing the observe kinematics to those expected as a function of the particle 
property of interest (e.g. m t ). The event probability density is estimated using a leading order matrix element and 
integrating over the unmeasured quantities. Detector effects are incorporated by integrating over resolution functions. 
The total likelihood for a given sample of events is obtained from the product of the individual event probability 
densities. Although it will not be discussed separately, it should be noted that the dynamic likelihood method is very 
similar to the matrix element method described here. Of all the methods explored at the Tevatron for determining m t , 
the matrix element method offers the best statistical sensitivity since it uses the full kinematic information available 
in each event and since it effectively gives more weight to events whose kinematics afford a more precise estimate of 
m\ ec ° by virtue of having a narrower event probability density. The principal downside is that it is enormously CPU 
intensive, typically requiring hours of CPU per event. This limitation renders these analyses less nimble than analyses 
using the template method. 

The method is based on an approach suggested in ( Dalitz and Goldstein . 1999t Kondo et all, 1993) and close to the 



meth od used for some of the FF-boson mass measurements at LEP ( Abreu et al . 1998 : Berends et al . 1998 : Juste) . 



1999). Currently, all the most precise m t determinations use the matrix element method. It has been applied in 



all the decay final states. In the Ijt and had channels it is extended to include an integration over the jet energy 
corrections, which are then constrained in situ using hadronic VF-boson de cays. A detailed des cription of the matrix 
element technique applied to measure the top-quark mass can be found in (|Fiedler et all |2010| ). 



The event probability P ev t is built from a tt probability (P t j) and a background probability (Pb), 

P GVt (x; mt,fu) = fit ' p tt( x '> m t) + fit) ' p ^( x ) ( 13 ) 

where x denotes the set of observed variables (i.e. the jet and lepton momenta) and fa is the it signal fraction in the 
event sample. The signal and background probability densities are constructed by integrating over the appropriate 
parton-level differential cross section, da(y)/dy, convoluted with parton distribution functions (k(q)) and resolution 
effects. The resolution effects are described by the transfer functions, W(x; y), which give the probability of observing 
the set of variables x given the underlying partonic quantities, y. Typically the jet and lepton angles are taken to 
be exactly measured so that the relevant transfer functions are simply Dirac delta functions. The transfer functions 
for the jet energies are parameterized as a function of parton energy and rapidity using fully simulated MC events. 
Separate transfer functions are derived for uds, gluon, and 6-jets. The transfer functions for the lepton energies are 
usually taken to be Dirac delta functions but are sometimes treated in a manner similar to the jets. It should be 
noted that the transfer functions are not assumed to be Gaussian, and more sophisticated functions can be used to 
obtain a more accurate description of the relevant resolutions. The it probability density is expressed as 

P ( X ) = ^L E / E ^-Hn)H<l2)dq 1 dq 2 W(x;y)dy (14) 

j — p comb flavors 

where the first sum is over all possible jet-parton combinations and the second sum is over all relevant PDF parton 
flavors. The probability density is normalized to the total observed cross section, <7 obs , after including event selection 
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effects. For signaf the differential cross section is taken from the leading order matrix element, Ai(qq —> tt — s- y), 
found in ( Mahlon and Parkel 1997 ). which depends on m t , — > p or background the differential cross 

section is taken from a sum of matrix elements, calculated by dedicated generators, and is independent of m t . 
The background probability density usually includes contributions from just the dominant background processes. The 
effect of the resulting approximation on the measured mt is included as part of the systematic uncertainties and is 
usually found to be quite small. For analyses that constrain the jet energy corrections using hadronic lU-boson decays, 
the jet energy transfer functions, and thus the event probability densities, are additionally expressed in terms of an 
overall jet energy scale factor, JES, W(x; y) — » W(x; y, JES). 

Both CDF and D0 have recently published m easurements of the top- quark mass in the Ijt channel using the matrix 
element technique. In their most recent analysis (jAaltonen et al . 20091 ) CDF also takes into account the finite angular 
resolution of the jets. In this analysis the gluon fusion diagrams are also included in the tt matrix element and the 
jet-parton combinations are weighted using information from a 6-jet identification algorithm. A NN discriminant is 
used to separate tt from background events rather than the matrix elements for the background processes. Events 
with very small likelihoods are vetoed in order to remove poorly behaved it events arising, for example, from the 
presence of ISR/FSR jets. The jet energy corrections are constrained using the W — > qq' decays. The method is 
calibrated using MC samples generated at various values of m t and assuming different JES factors. Using 1.9 fb _1 
of data CDF measures m t = 172.7 ± 1.8 (stat + JES) ± 1.2 (syst) GeV/c 2 (see Fig. IT5) . The dominant systematic 
uncertainty arises fr om difference between MC generators. 

The D0 analysis ( Abazov et all . l2008eh includes only the qq — > tt component in the signal matrix element and 
calculates the background probability density using only the matrix element for the W + 4 partons process. The jet- 
parton combinations are weighted using information from a 6-jet identification algorithm and an in-situ JES calibration 
is performed. The method is calibrated using MC samples generated at various values of m t and assuming different 
JES values. Using 1 ftT 1 of data D0 measures m t = 171.5 ± 1.8 (stat + JES) ± 1.1 (syst) GeV/c 2 (see Fig. [15). The 
dominant systematic uncertainty arises from uncertainties associated with calibrating the response of the calorimeter 
to 6-quark jets. 



CDF has also published a measurement using the matrix element technique in the dil channel (jAaltonen et 



2009j) . In this channel, no JES in-situ calibration is possible. In thi s analysis, the event select i on has been optimized 



to get the best statistical precision on m t using neuroevolution ([Stanley and Miikkulainenl, l2002h to train a NN 
discriminant. The signal matrix element includes only the qq — ¥ tt component. Background probability densities are 
computed using the matrix elements for the Z/j ee, /i^+jets, W+ > 3 jets where a jet is misidentified as a lepton, 
and IFIU+jets processes. Events with or without identified 6-jets are fit separately. Using 2 ftT 1 , CDF measures 
mt = 171.2 ± 2.7 (stat) ± 2.9 (syst) GeV/c 2 . The dominant systematic uncertainty is due to uncertainties in the jet 
energy corrections. 



Ideogram Method 

The ideogram method can be thought of as an approximation to the matrix element method. A kinematic fit is 
used to determine m\ ec ° for each event. A per event probability density is calculated as a function of m t by calculating 
the probability of observing m\ cc ° assuming the true top-quark mass is mt and knowing the resolution o m ™ ca ■ All 
jet-parton combinations are considered for each event and each carries a weight derived from its corresponding fit x 2 ■ 
Typically, additional weights are included that use information from 6-jet identification algorithms. Like the matrix 
element method the event probability density is built from a signal and a background piece. The signal probability 
density is a convolution of a Gaussian with a Breit-Wigner while the background probability density is taken from 
MC simulation. The total likelihood for a given sample of events is obtained from the product of the individual event 
probability densities. A JES constraint can be incorporated by repeating the kinematic fits for different assumptions 
on the jet energy corrections scale factor. Like the matrix element method, the ideogram method offers improved 
statistical sensitivity. In terms of sophistication, CPU budget, and statistical sensitivity the ideogram method typically 
falls between the template and matrix element methods. It should be noted that since it requires a full reconstruction 
of the tt kinematics, this method can only be used in the Ijt and had final states. 



D0 has performed a measurement using this technique in the Ijt channel ([Abazov et all [2007cj). The event 
probability density is factorized into the product of two separate probabilities for signal and background each. The first 
probability depends on m t and JES and provides the necessary information to constrain the jet energy corrections and 
to measure the top-quark mass. The second probability depends on the output of a multivariate discriminant designed 
to separate tt from background processes. The variables used in the discriminant are chosen and the discriminant 
itself is constructed to be uncorrelated with m t and JES. This second probability helps constrain the observed tt 
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fraction in the event sample and to de-weight events that have kinematics more consistent with having originated 
from background processes. MC simulation is used to account for contributions from wrong jet-parton assignments 
and background events. Using 0.43 fb -1 , D0 measures m t = 173.7 ± 4.4 (stat + JES)t^ (syst) GeV/c 2 . The 
dominant systematic uncertainty arises from uncertainties associated with calibrating the response of the calorimeter 
to 6-quark jets and from uncertainties associated with modeling the tt signal. 

The method was used by CDF to produced a result in the had channel using 310 pb" 1 (|Aaltonen et all . 12007). The 
result was limited by the uncertainties in the jet energy corrections. As described above subsequent m t determinations 
in the had channel exploit W — ¥ qq' decays to reduce the jet energy correction systematic and also include significantly 
larger data sets. 



Using the tt Cross Section 



As mentioned above, there is some ambiguity associated with the theoretical interpretation of the mt parameter 
measured by the above techniques. Assuming SM production and decay, an estimate of the top-quark mass can also 
be made by comparing the measured production cross section to theory predictions. While less precise, this method 
has the advantage that the m t parameter in the predictions is theoretically well defined. 

D0 has performed this measurement ([Abazov et all l2009q) . The signal acceptance is estimated as a function of 



mt using MC samples generated at various m t values. The resulting acceptances are smoothed and the production 
cross section is then measured as a function of mt as shown in Fig. [13j The m t dependence of the acceptance arises 
from the ev ent selection lepton and jet Et requirements. The measured cross section is compared to several NLO 
predictions (jBeenakker et al , Il989t ICacciari etall . 120081: iMoch and Uwerl . 120081 iNadolskv eTaR 120081) that define m t 



as the pole mass. A normalized likelihood function is formed as a function of mt by comparing the theory prediction 
with the measured cross section at each assumed top-quark mass. The uncertainty in the theory arising from variations 
in PDFs and in the renormalization and factorization sca l es an d the uncertainty in the measured cross section are 
accounted for in the comparison. Using (IMoch and Uwerl . l2008h and the combined Ijt + dil measured cross section 
from 1 fb~ of data, D0 measures mt — I69.1I52 GeV/c 2 . This result is in agreement with the direct top-quark 
mass measurements. Using the same theoretical predictions but an expected experimental uncertainty on the it cross 
section of 6% at the end of the Tevatron run, a top-quark mass uncertainty of around 3.3 GeV/c 2 could be achieved 
with this method. It should be n oted that an extraction of the running mass in the MS scheme have also been 
performed ( Langenfeld et al . 20091 ). 



Combination 



The world average top-quark mass is obtained from a combination of CDF and D0 most precise mt measurements 
using up to 5.6 pb _1 of data. Five measurements using data taken in an earlier run (1990-1995) and six measurements 
using data from the ongoing run (2001-) are used. Work has been performed between the two collaborations to stan- 
dardize the assessment of systematic uncertainties. For the combination a detailed breakdown of the various sources of 
uncertainty has been established in order to properly account f or correlations among the various measurements. The 



19881). The prelimina ry combined value for the top-quark 
173.32 ± 0.56 (stat) ± 0.89 (syst) GeV/c 2 (|Aaltonen et all , l2010al ) which corresponds to a total uncer- 



combination uses the method described in (jLvons et 

mass is to* 

tainty of 1.06 GeV/c 2 and a relative precision of 0.61 %. The combination has a x 2 of 6.1 for 10 degrees-of-freedom 
indicating good consistency among the measurements across different experiments and different decay channels. The 
most recent CDF and D0 measurements in the Ijt channel using the matrix element method carry the largest weights 
in the combination followed by the CDF template measurement in the had channel. 

In the future, the statistical and in-situ JES related uncertainties will continue to decrease with the statistics of 
the data samples employed. The systematic uncertainty has roughly equal contributions from residual jet energy 
correction uncertainties and uncertainties associated with modeling the tt signal from variations in MC generators 
and color reconnection effects. Work is in progress to better understand these sources of systematic uncertainty with 
the aim of reducing them. A total top-quark mass uncertainty below 1 GeV/c 2 is potentially achievable in the near 
future. 



3G 



2. Mass Difference 



The triple product of the discrete symmetries of charge conjugation (C), parity conjugation ( P), and time re- 
versal (T) is an exact symmetry of any local Lorentz-invariant quantum field theory and requires that the par- 
ticle iajid_antijDart^ However, CPT can be violated in some non-local new physics mod- 
els (IColladay a nd Kostclcclo; ,_ 19971 1 . Tests of CPT have been performed in many sectors with no evidence of CPT 
violation ( Nakamura et aZ.I . l201dh . Measuring a q — q mass difference is usually difficult since hadronization introduces 
QCD binding and non-perturbative showering effects. However, since it decays before hadroniz ing the top quark 
provi des a unique opportunity to measure a q — q mass difference free from these QCD effects (jCembranos et al. 
20081) . 

D0 performed the measurement of the t — t mass difference in t he Ijt final sta t e using 1 fb -1 of data (lAbazov et al. 



2009dl ). This analysis uses the matrix element technique as in (jAbazov et all l2008ef) to form probability densities 



for each event as a function of the mass of the t and the t and the tt signal fraction. SM it production and decay 
is assumed and the lepton charge is used to determine whether the ivb invariant mass should be assigned as rat or 
raj. Many of the systematic uncertainties relev ant for the top-quark mass measurements cancel here. The result is 

nit — "mi = 

3.8 ± 3.4 (stat) ± 1.2 (syst) GeV/c 2 (jAbazov et a/.l . l2009dh . consistent with CPT invariance. 

-l 



A new CDF measurement using the template method in the Ijt channe l with 5.6 fb 
-3.3 ± 1.4 (stat) ± 1.0 (syst) GeV/c 2 , also consistent with CPT invariance (jAaltonen et 



of data finds rat — TOj 
2011dl) . 



3. Charge 

The electric charge is a fundamental property of the top quark predicted in the SM to be +2/3 in units of the 
electric charge of the electron. An experimental determination of the top-quark electric charge serves as a check 
of SM consistency. In all the tt analyses so far discussed the assumption is made that t —> W + b and t —> W~b. 
However, since none of these analyses have attempted to differentiate a 6-jet from a b — jet, it is possible that what 
is being observed is an exotic top-like particle that decays to W + b and W~b. Just such a particle can be found 
in (jChang et al. , Il999t IChoudhurv et all l2002jb a four generation model that postulates a —4/3 charged quark that 



decays to W~b. In this model, the right-handed 6-quark mixes with a —1/3 charged heavy quark Q\. Its weak isospin 
partner Q4 has charge -4/3, a mass around 172.5 GeV/c 2 and would be the particle observed at the Tevatron while 
the left-handed real top quark would have a mass around 270 GeV/c 2 and would thus far have escaped detection. It 
is worth noting that the precision electroweak data can be satisfactorily described by this model. 

The electric charge of the top quark can be determined from measurements of the electromagnetic coupling strength 
using ttj events or by reconstructing the electric charges of the top-quark decay products. The production cross section 
for pp — > ttj is proportional to the square of the top-quark electric charge. However this measurement suffers from a 
large irreducible background coming from processes where the photon is radiated from the 6-quark or the W boson. 
Both CDF and D0 use the second approach to measure the top-quark electric charge in the Ijt channel. A kinematic 
fit is used to fully reconstruct the tt kinematics and choose particular jet-parton and 6-jet pairings (ie. is the 6-quark 
jet paired with the leptonically decaying W, ivb, or the hadronic decaying W, qq'b). The flavor of the 6-jets must 
then be specified by determining whether the jet originated from a 6-quark or a 6-quark. The sum of the electric 
charges of the charged lepton and the appropriate electric charge of the 6-jet for the Ivb — jet pairing is used as an 
estimate of the top-quark charge. 

CDF requires that at least one jet is identified as a 6-jet using an algorithm that exploits the long lifetime of 
6-hadrons. A kinematic fit is employed and the jet-parton and 6-jet pairings are taken from the permutation yielding 
the lowest fit chi-squared. This corresponds to the correct assignment about 76% of the time. Soft lepton taggers 
identify the flavor of the 6-jets with a 69% purity. The purity of the 6-jet flavor determinati on is calibrated using 66 



data and MC events. Using 2.7 fb - the CDF data is consistent with the SM expectation (jAaltonen et all l2010bj ) 
and excludes the exotic -4/3 charged quark interpretation at 95% confidence level. 

D0 requires that at least two jets are identified as 6-jets and also uses a kinematic fit to reconstruct the tt kinematics 
and choose the jet-parton and 6-jet pairings. Using the permutation with the lowest fit chi-squared yields the correct 
assignments about 84% of the time. A jet-charge algorithm is used to identify the flavor of the 6-jet. The optimized 
estimator of the jet charge is qj et = QiPTi' 6 ) / (J2iPTi' e ) where i runs over all tracks associated with the jet in 
question. The expected qj e t distributions for 6, 6, c and c jets are derived from di-jet data samples enhanced in heavy 
flavor. The charge of the leptonic decaying top quark is reconstructed as the sum of lepton charge and the charge of 
the 6-jet assigned to this top quark while the charge of the hadronic decaying top quark is computed as the sum of the 



second 6-jet minus the lepton charge. Using 0.4 fb 1 the D0 data is consistent with the SM expectation ( Abazov et all 
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2007bl ) and excludes the exotic -4/3 charged quark interpretation at the 92% confidence level. 



4. Width 



The total width of the top quark is precisely determined in the SM once m t is specified. Because m t is large, so 
is the top-quark width, r t , and consequently its lifetime, ~ l/r t , is extremely short, which precludes measuring the 
width using reconstructed displaced vertices from t — ¥ Wb decays. A measurement of T t can also be determined from 



the width of the m\ ec ° di stribution. In the SM at N LO, neglecting terms of order m 2 /m 2 , a s , and or s /i 



top-quark total width is (jJezabek and Kuhn . Il989h : 



/m 



I the 



n i 



HV 



my 



2a s /2tt 2 
"3tT l~3~ 



(15) 



where T® 



G F my] 
87rV2 



Vtb\ 2 with Gf being the Fermi coupling constant. This yields approximately to T t — 1.4 GeV/c 2 

for a 172 GeV/c 2 top-quark mass. Deviations from the SM expectation could signal the precense of BSM physics 
contributions arising, for example, from t — > H + b decays, anomalous tWb form factors, enhanced rates of t — > dW + 
or t — >• sW + decays, or the existence of a fourth-generation fo'-quark. 

CDF uses 1 fb _1 of data in the Ijt channel with at least one identified 6-jets and employs a kinematic fit to 
obtain the m\ cco distribution. The shape of the m\ eco distribution is sensitive to the top-quark width and a tem- 
plate meth od is employed to measure T t . The CDF data is shown in Fig. [14] A confidence interval based on the 
method of ( Feldman and Cousinsl 1998) is used to determine T t < 13.1 GeV/c 2 at 95% confidence level. The dom- 
inant syste matic uncertaint i es aris e from uncertainties in the jet energy corrections and in modeling the jet energy 



resolutions ( Aaltonen et al . 2009b) 



j eco it is difficult to reach 



Since the SM top-quark width is far smaller than the experimental resolution on m 
high precision using the direct measurement technique described above. However it is also possible to extract T t by 
determining its partial width T(t — > Wb) using the measured cross section for the electroweak pro duction of single top 
quarks and the measure d top-quark branching fraction B(t — > Wb), T t — T(t —> Wb)/B(t — > Wb) (jCarlson and Yuan . 
19951: lYuanl . Il993l Il995h . This method assumes that the couplings at the tWb production and decay vertices are the 
same, that the production via flavor changing neutral currents is negligible, and that |Vtd|, \Vta\ « \Vtb\- The 
B(t — ► Wb) is determined from the Rb measurement discussed in Sec. IIX.AI assuming B(t — > Wq) = 1. Electroweak 
single top production can occur via an s-channel or a t-channel exchange of a W boson. If the s + i-channel cross 
section is used to extract r t , then the SM ratio (s/t) is assumed. To minimize the model dependence of this method, 
the dominant ^-channel cross section can be used instead. 

D0 uses the 1 fb _1 Rb measurement (Abazov et all 2008h) and a 2.3 fb _1 t-channel single top-quark production 
cross section measurement ( Abazov et aLl . l2010el) to determine the top-quark width using this method. A NLO QCD 
calculation is used to determine B(t —> Wb) from the t-channel cross section and the correla t ed syst ematic uncertainties 



have been taken into account. Using these inputs the top-quark width (jAbazov et all l2011d) is measured to be 
r t = I.99+0 55 GeV/c 2 in agreement with the SM expectation (cf. Fig. [T4")) . This top-quark width measurement can 
be used to extract a limit on | Vtv \ , an element of a four generation CKM matrix characterizing the coupling strength 
of the top quark with a hypothetical fourth-generation 6'-quark. Assuming my > mt — m w an d the unitarity of the 
extended CKM matrix, this measurement yields \Vtb' \ < 0.63 at 95% confidence level. 



5. Spin Correlations 



At the Tevatron the spins of the t and t at production are correlated (jBarger et al. . 119891: iMahlon and Parkelll996h . 
The top quarks decay on a time scale that is short relative to the hadronization and spin-flip time sca le and, hence, the 
spin correlations are preserved and affect the angular distributions of the top-quark decay products ( Falk and Peskinl . 
1994J). If we note 9 the angle between the spin quantization axis and the direction of flight of the VF-boson decay 



particle in the rest frame of its parent top quark, calculations show that the angular co rrelations are strongest fo r 
the charged leptons (9i) or the down-type quarks (9d) from the decay of the W bosons ( Mahlon and Parka l2010h . 



Experimentally, the two dimensional distribution (cos6^, cos 9d) or (cos (9^+, cos 6^-) is used to determine the spin- 
correlation coefficient, C, 
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where a denotes the total tt cross section in the relevant decay channel, and the indices 1 and 2 re fer to the top-quark 
and a ntitop-quark decay products respectively. This relationship has been calculated to order (jBernreuther et < 



2004al ). The angles are measured relative to a given axis. Three common choices are a) the beam basis, which uses the 



direction of one of the colliding hadrons as the reference axis, b) the helicity basis, which uses the t flight direction as 
the reference axis, and c) the off-diagonal basis, which uses a ref erence axis chosen so that the contribution from like- 
spin t — i pairs vanishes at leading order in perturbation theory (jParke and Shadmil . If 9961 ). The measured coefficient 
changes as a function of m t j and depends on the fraction of tt pairs produced at threshold and the fraction produced 
by gluo n-fusion. In the absence o f selection or acceptance requirements C can be computed in terms of expectation 
values ([Bernreuther et all l2004al) C = — 9 < cos#i cos 62 >, which vary with the choice of basis employed as shown 



in Table HU 

The observation of t — i spin correlations would confirm that the top quark is a spin 1/2 particle and that it 
decays before strong interaction processes affect the t and i spins. Several BSM mo dels could affect the measured 
coefficient. The presence of a techni-77 in two scale technicolor ( Eichten and Lane]. 1994 ) would enhance the unlike spin 
contributions. A new vector particle associated with topcolor ( Hill and Parkei 19941) would affect the spin correlations 
by changing the fraction of gluon-fusion produced it pairs. On the decay side, BSM contributions from t — > H + b 
decays, for example, would affect the cor relations between the top- quark spin and its decay products and therefore 
modify the observed angular correlations (jMahlon and Parkei Il996l) . 

Both CDF and DO have performed preliminary measurements of the spin correlation in the dil channel using 
2.8 fb _1 and 5.4 fb~ of data, respectively. To fully reconstruct the tt kinematics CDF uses a likelihood method 
to determine the best solution for the unknown neutrino momenta using the observed lepton, jet, and S^t energies 
and taking into account their relevant resolutions, while D0 uses the neutrino weighting technique described above 
in Sec. IIX.C.11 After reconstructing the tt final state, two dimensional (cos 6^-, cos6^+) templates are constructed. 
The signal templates are built by reweighing Pythia tt samples using Eq. Q1)J A maximum likelihood fit to the 
data, including systematic uncertainties via nuisance parameters, gives C = 0.32^°}^ f or CDF in the off-diagonal 
basis (jAaltonen et all l2009~fl ) and C = 0. 10+^45 for D in the beam basis (|Abazov aZ.L l2011ah . 

CDF has also performed a measurement of the tt spin correlation in the helicity basis using the Ijt channel in 
4.3 fb _1 of data. The difficulty here lies in determining which jet comes from the hadronization of the down-type 
quark. It is assumed to be the jet cl osest to the &-jet in the_W-boson rest frame. This prescription gives the right 
assignment about 60 % of the time ( Mahlon and Parkei Il996l ). A templ ate fit to the two dime nsion distribution 
(cos0£, cosOd) is used to determine C = 0.60 ± 0.50 (stat) ± 0.16 (syst) (jAaltonen et ai , 2011bl) . Using the same 
methodology CDF updated the analysis to use both the helicity and beam bases and to include 5.3 fb _1 of data. The 
preliminary results are C = 0.48 ± 48 (sta t) ±0.22 (syst) in the helicity basis and C = 0.72 ±0.64 (stat) ±0.26 (syst) 
in the beam basis ( Aaltonen et all l2010dr ). One dimensional distributions of cosQi ■ cos, Od.i from the DO and CDF 
analyses are shown in Fig. 1151 

All these results are completely dominated by the statistical uncertainty and will benefit from including the full 
Tevatron statistics in the future. This measurement will remain interesting at the Tevatron since a different spin 
correlation will be produced at the LHC where tt production is dominated by gluon-gluon fusion rather than qq 
annihilation. 



X. DIRECT SEARCHES FOR NEW PHYSICS IN THE TOP-QUARK SECTOR 

A. Search for tt Resonances 



Many SM extensions predict massive new p articles that can d ecay to a tt pair. For example, the scalar color- 
singlets predicted in two- Higgs double ts models ( Dicus et ai , 19941) . the vector color-singlet Z' par ticles postulated i n 
extended gauge theories ( LeikeL 19991 ). the KK g \ uoa and KKz states of extra dimension mode ls (Lillie et a/J. 20071). 
and the axigluons, colorons, and vector color-octet particles that appear in top-color models ( Hill and Parkelll994t 



Sehgal and Wanningerl . 119881 ). all couple to tt pairs and would affect the m t j distribution. In particular, resonances 



whose natural width is small relative to the experimental resolution would appear as an anomalous peak. Both CDF 
and DO have searched for such peaks using the Ijt channel. In all such searches so far performed at the Tevatron, the 
natural width of the resonance, generically denoted Z' since Z' models are often used as benchmark models, is taken 
to be Tz> = 0.012mz', well below the typical experimental resolution achieved of 0.05 — 0.10 x m t j. The resulting 
upper limits derived from these analyses are, in general, valid for all widths up to Y z> ~ 0.05mz> unless otherwise 

specified. 

Different approaches were used to reconstruct m t j for each event. CDF uses (jAaltonen et all l2008hl) a matrix 
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element technique similar to the one developed for the top-quark mass measurement as described in Sec. IIX.C.11 The 
mean of the event probability density distribution is used to determine the reconstructed m t j value for that particular 
event. Templates are then formed for a variety of postulated Z' resonances and including background contributions 
from ti, W+jets, and QCD processes. A binned likelihood fit to the data is performed letting the Z' production cross 
section vary in the fit. Using 0.68 fb~ , the data is found to be consistent with SM expectations and C DF excludes 



a leptophobic top-color resonance candidate with a mass less than 725 GeV/c 2 (jAaltonen et all l2008h|) at the 95% 
confidence level. A complementary analysis uses a slightly different event selection and a kinematic fit to obtain 
the reconstructed m t j value for ea ch event. With a 0.96 fb -1 data set this CDF analysis excludes a leptophobic Z' 
with a mass less than 720 GeV/c 2 ( Aaltonen et al , 2008el) at the 95% confidence level. The production cross section 
for any narrow Z'-like new particle decaying to ti is determined to be less that 0.64 pb at the 95% confidence level 
for all mz> > 600 GeV/c 2 . All the above limits include the effects of systematic uncertainties on the background 
normalizations, the signal efficiencies, and the shapes of the signal and background m t j distributions. 

D0 has performed a similar search in the Ijt channel requiring at least one identified fo-j et in each event. A kinematic 
fit is employed to determine the reconstructed m t j value for each event and a fit to the data is performed using a binned 
likelihood based on the template method. No significant deviation from the SM expectation is observed, as shown in 
Fig. [Tni an d a Bayesian approach is used to set 95 % confidence level limits on ax • B(X — > it) for different values 
of mx ■ Systematic uncertainties that affect both the normalization and the shape of the signal and background m t j 
distributions are taken in to account. Using 0.9 fb _1 D0 excludes a lepto-phobic Z' with a mass less than 700 GeV/c 2 



at 95 % confidence level (jAbazov et ali\2008m . 



CDF has also searc hed for a massive vecto r color-octet particle, G, such as a massive gluon that preferentially 



couples to top quarks (jAaltonen et all l2010il ). Neglecting the gluon- G coupling, the strength of the massive gluon 
coupling to any quark is X q a Sl where a s is the SM strong coupling constant. In this analysis the mass and width 
of the particle G are denoted M and T. If only SM top-quark decays are assumed, the ti production through an 
intermediate G interferes with the SM ti production. After selecting Ijt events the m t j spectrum is reconstructed 
using a dynamical likelihood method similar to the matrix element methods described in Sec. IIX.C.11 A standard 
unbinncd likelihood technique is then applied to fit the data. If the mass M and T/M are fixed, the signal m t j 
distribution depends only on the assumed coupling strength At. Taking into account systematic effects, limits on At 
as a function of M for T/M = 0.1 and T/M = 0.5 are determined. Using 1.9 fb _1 the data, shown in Fig. [T6l are 
consistent with SM expectations and CDF excludes the presence of massive gluons with couplings to the top quark 
larger than 0.5 for 400 < M < 800 GeV/c 2 and 0.05 < T/M < 0.5 at the 95% confidence level. 



B. Search for Top-quark like t'-Quarks 



A number of BSM models predict a fourth generation of massive fermions (jFrampton et all 2000; iHoldom et al 
2009). The precise measurement of the Z-boson invisible width excludes the existence of a fourth neutrino family with 



mass less than half the Z-boson mass but a new neutral lepton with a mass larger than that is still comp atible with 
electroweak da ta. An extra ch iral fermion family appears in gr and-unified theories (GUTs) (|Rossl . Il984f ) or TV = 2 



SUSY models (iHe et al 
partner b' is preferred so that for nit 



20011 ). Often a small mass difference between the fourth generation t' and its down- type 
> mi,' + mw and moderate mixing between t' and the SM light down-type 



quarks q, the decay mode t 1 — > Wq dominates. 

In case the fourth generation fermions are vector-lik e, experimental constraints are even weake r. Such heavy exotic 



20031 ) or models with large 



vector-like quarks are predicted by little Higgs models (jDobrescu and Hiil . [l997 
extra-dimensions. If not too heavy, such particles can be pair produced at the Tevatron and can decay to Wq with 
a substantial branchi ng fraction. Some mod els predict the production of a single vector-like heavy t' in association 



2009) 



with a SM top quark ( Dobrescu et all 

Using 0.76 fb~ , CDF has searched f or the pair production of massive up-type t! quark/anti-quark pairs that 



each decay to Wq using the Ijt channel (jAaltonen et all l2008gf ). The analysis employs a kinematic fit to determine 
TO reco f or eacn even t anc j uses the distribution to further discriminate a possible signal from SM processes. A 
two dimensional binned likelihood fit of the background and signal to the observed (Ht , m T t f c °) distribution in data 
is performed to set limits on the t' mass. Systematic uncertainties are treated as nuisance parameters and include 
variations in signal and background normalizations and shape. While still statistically limited systematic uncertainties 
are dominated by the uncertainty on the jet energy calibration. In the fit the ti production cross section is constrained 
to its SM value. Assuming B(t' — > Wq) = 1, CDF excludes a t' mass below 256 GeV/c 2 at 95 % confidence level. 

CDF and D0 have provided preliminary updates for this t' search with more luminosity using similar analysis 
techniques. Using 5.6 fb~ , CDF excludes t' with masses below 358 GeV/c 2 , while D0 using 4.3 fb _1 excludes t' 
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masses below 296 GeV/c 2 both at 95 % confidence level ([Aaltonen et all 1201 let lAbazov et a^ , |2010a) 



C. Search for Scalar Top Quarks 



In supersymmetric models ( Haber and Kane , 1985 ; Martin , 1997 ; Nillesl . 1984 ) each of the SM fermions has a scalar 
partner. In particular the left- and right-handed top quarks each have a scalar super-partner: tn and t^. The mass 
eigenstates t\ and ti of these two scalar top quarks are rotated with respect to tn and t^. Due to the large top-quark 
mass, this mixing can be substantial and t\ and ii can have a large mass difference. Hence, the lightest scalar top 
quark can be the lightest scalar quark and even lighter than the top quark itself. In the MSSM the scalar top quarks 
are mainly produced in t\t\ pairs via the strong interaction. In general the production cross section is relatively small. 
For example, when mi — rat the t\t\ pair production cross section is about 10 times smaller than the it production 
cross section. The t decay differs according to the parameters of the model and is particularly dependent on the 
masses of the super-particles involved. If kinematically possible, it can decay via the two body processes t — > txi, 
t — s- bxl , or i — > cxi, where Xi is the lightest neutralino and xt the lightest chargino. Three-body decays are also 
possible, t — > W + bxi, i — > b£ + vi, or t — > b£ + vi, where £ and i>i are the super-partners of a given lepton I and neutrino 
vg. For scalar top quarks lighter than the top quark the decay channel t —> xt^ — > (XiW + )b can dominate. In 
this decay channel the titi pair signature is similar to the ti one apart from the presence of the neutralinos, whose 
experimental signature mimics that of a neutrino and are mostly produced back-to-back. 

D0 searched for some admixture of t\ti pairs in the ti final state using the Ijt channel (jAbazov et 



this analysis the neutralino mass has been fixed to a value close to the experimental LEP limit, m^o 



I2009gj) . In 

50 GeV/c 2 , 



while the t mass has been varied from 130 to 190 GeV/c 2 and the chargino mass from 90 to 150 GeV/c 2 . Events with 
at least three jets and at least one identified fo-j et have been used. For events with at least four jets, a kinematic 
reconstruction to the it hypothesis is performed. As the it signature is very similar to the t signal, a multivariate 
likelihood discriminant is employed to separate the two processes. Separate discriminants have been trained for the 
three-jet and the four-jet events. Variables based on the two jet or three jet invariant masses, the reconstructed top- 
quark mass, the jet Et, the scalar sum of the jet Et, and the angular separation between the jets are used as inputs 
to the discriminants. As no excess over the SM expectation is observed, a binned likelihood over the discriminant 
distributions is constructed using Bayesian statistics to extract the upper limits on the t\t\ pair production cross 
section taking systematic uncertainties into account. Using 0.9 fb _1 , D0 obtains limits that are a factor of 2 to 13 



times larger than the MSSM theory prediction and agree with the expected limits within uncertainties (jAbazov et al. 
20093). 

When both scalar top quarks decay as t — ¥ blvi the experim ental signature is simila r to the it dil signature. Searche s 



exploiting this decay channels have been performed by CDF ( Aaltonen et al. . 2010ll) and D0 ( Abazov et aZ.I . [2008g^) 
A similar experimental signature is obtained if the scalar top quark is assumed to decay as t 



Xib ~ ► Xi£ +l/ b as 



is done in ( Aaltonen et al . 2010kt ). In all three of these analyses the selection criteria are optimized to suppress the 



tt contributions. The data are consistent with the SM expectation and the resulting limits significantly extend the 
LEP 2 exclusions when >> m^m^o. 



XI. PROSPECTS AND CONCLUSIONS 

Owing to its large mass, the top quark was discovered relatively recently and potentially offers a unique window 
into Beyond-the-Standard-Model physics. Using data samples with C(10 3 ) candidate ti events, the CDF and D0 
experiments at the Fermilab Tevatron have pioneered numerous analyses to exploit these unique events. They have 
instituted a broad physics program that tests the Standard Model descriptions of the ti production mechanisms, the 
top-quark decay widths, and the intrinsic properties of the top quark such as its electric charge and mass. They have 
also probed the it event samples for direct evidence of new particles that might mimic the experimental signature of the 
top quark, such as the scalar top quarks of Supersymmetry or the fourth generation t'-quark of Grand Unified Theories. 
A few of the measurements are systematics limited and have reached precisions comparable to those associated with 
the relevant theory predictions, such as the measured production cross section and the measured top-quark mass. 
The data provide strong evidence that the top quark carries electric charge +2/3 in units of the electron electric 
charge. Measurements of the polarization fractions of W bosons from top-quark decay and of it spin correlations 
are consistent with the Standard Model expectation of a spin 1/2 particle with a V — A coupling to the charged 
weak current. The decay widths of the top quark are also consistent with the Standard Model expectations, being 
dominated by t — ¥ W + b and with no evidence of anomalous flavor-changing-neutral-current or charged higgs decays. 
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The results of these many measurements when compared across the various ti final states yield consistent results - 
again as expected by the Standard Model. In general, the detailed exploration of the it sample so far performed 
by the Tevatron experiments reveals no significant discrepancies with the Standard Model. Except for a few cases, 
most all the measurements are limited by their statistical uncertainty. So the precise study of the top quark will 
persevere as the Tevatron is continuing to deliver collision data. The LHC experiments expect to collect it samples 
that are a few orders of magnitude larger in the coming several years. With these sorts of statistics most all the 
measurements should achieve levels of precision that more thoroughly probe the Standard Model expectations and 
more stringently restrict the parameter space of a wide variety of new physics models. Even for those measurements 
already systematics limited at the Tevatron, the LHC experiments may be able to exploit the larger statistics to 
identify a sub-set of events less affected by the systematics in question or to employ analysis techniques which trade 
statistical uncertainties for systematic uncertainties to achieve an improvement in the total uncertainty. The LHC is 
poised to build upon the foundation laid by the Tevatron experiments, and use the physics of the top quark to more 
fully elucidate the workings of nature at its most fundamental level. 
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FIGURES 



FIG. 1 Example of Feynman diagrams depicting SM interactions for the strong (left), electromagnetic (middle), and weak 
(right) forces. 



FIG. 2 Feynman diagrams of pp — > tt for the qq annihilation process (left) and the gluon-gluon fusion process (middle and 
right). 




FIG. 3 Feynman diagrams of pp — > tt production with the tt decaying to the dilepton (left), lepton+jets (middle), and all 
hadronic (right) final states. 
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FIG. 4 A cut away view of the CDF (top) and D0 (bottom) detectors. 



47 





D0, L=5.3 lb'' 




1 >2 
Number of b-tagged jets 



FIG. 6 The predicted and observed number of events as a function of jet multiplicity for CDF (left) and as a function of 
the number of identified fa-jets for D0 (right) bo th from a it cross section analysis in the Ijt channel using &-jet identifica- 
tion (|Aaltonen et all l2010cl : lAbazov et oil l2011dl ') . 
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FIG. 7 The observed S^t distribution compared to the SM e xpectations for the tt c r oss section measurem ent using CDF dil 
events (left) and using D0 11 and It events combined (right) (|Aaltonen et al. l l2009al : lAb azov et aUl2009j ). 
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FIG. 8 The NN output distribution for QCD data (solid) and ti MC events (dashed) from CDF (left). A comparison of 
the likelihood distribution between data and MC using the signal purity fit fr om the data from D0 (right). Fits to thes e 
distributions are used to measure the ti cross section in the all hadronic channel l|Aaltonen et all l2010el ; lAbazov et all I2010J ) . 
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FIG. 9 The measured do/dm t j from CDF (left) compare d to the SM expectation. The measu r ed da /dpr from D0 (right) 
compared to QCD calculations and various MC generators (jAaltonen et all l2009d : lAbazov et al 
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FIG. 10 The CDF determined Aft in the tt frame, corrected to the parton level, compared to NLO predictions in bins of the 
ti invariant mass (left). The 95% confidence level limits on the fraction of ti produced via a Z' resonance as a function of the 

Z' mass (right) obtained from the D0 A 1 ^" 1 " 1 measurement (jAaltonen et al I l2011al ; lAbazov et all l2008d\ 
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FIG. 11 The rrit cco distribution (left) in the CDF Ijt channel with the signal and backgr ound template estimat es overlaid using 
rat = 172 GeV/c 2 and the nominal jet energy corrections for events with 1 identified 6-jet (jAaltonen et a l.. 2009c) . The likelihoo d 
curve from the data fit (right) in the Ijt channel using Lt and lepton Et as top-quark mass estimators l| Aaltonen et ql.U2010hl ). 
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FIG. 12 Results from mt measurements in the Ijt channel using the matrix el ement method. The 2-D likelihood (left) fro m 
CDF data and the data likelihood projected onto the m t axis (right) from D0 (jAaltonen et al 1, 1200911 : IXb azov et aZ.L l2008eh . 
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FIG. 13 D0 measu rement of the top-qua rk mass determined by comparing the measured and the SM predicted tt cross section 
as a function of m t (|Ab azov et oU l2009d ). 
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FIG. 14 The m T t cc ° distribution (left) for CDF events with identified 6-jets overlaid with signal and b ackground distributions 
and the poste r ior pro bability density (right) for the expected and measured Tt determination from D0 (jAaltonen et al 2009b; 
lAbazov et aILl2011ch . 
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FIG. 15 The CDF cos #4 • cos #4 distribution (left) in Ijt events using the beam basis. The D0 cos6 £ - ■ cos6 £ + distribution 
(right) in dil events using the beam basis. 
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FIG. 16 A comparison of t he observed mg distribution (points) to th at expected from SM sources (stacked histograms) for 
CDF (left) and D0 (right) l|Aaltonen et all l2010il ; lAbazov et aUl2008tf ). 



TABLES 



TABLE I The measured masses of the elementary particles in the Standard Model. For quarks other than the t op, it is estimated 
from hadron masses. For the neutrino masses only limits exist. Those listed are at the 90% confidence level (|Nakamura et all 
l201Ch . 



quark masses 


m u 


~ 0.002 


m c = 1.27 m t = 171.2 


(GeV/c 2 ) 


ma 


~ 0.005 


m s = 0.104 m b = 4.20 


(GeV/c 2 ) 


lepton masses 




< 2 x 10 


" 6 eV/c 2 , Yji m li > 2 x 10 


- 17 (eV/c 2 ) 2 


m e 


= 0.511 


m M = 105.6 m T = 1776.8 


(MeV/c 2 ) 



TABLE II NLO computation for the spin-correlation coefficient, C, discussed in Sec. IIX.C.5I ( Bernreuther et aH [2004a) 



basis 


dilepton channel 


lepton+jets channel 


helicity 


-0.352 


-0.168 


beam 


0.777 


0.370 


off-diagonal 


0.782 


0.372 



